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EXECUTIVE SUMMARY 

This Phase II (1994) Annual Progress Report presents two major report - 
sections describing the . thermal analysis of tilting- and flexure-pad hybrid 
bearings, and the unsteady flow and transient response of a point mass rotor 
supported on fluid film bearings. A literature review on the subject of 
two-phase flow in fluid film bearings and part of the proposed work for 1995 are 
also included. 

The programs delivered at the end of 1994 are named hydroflext and 
hydrotran. Both codes are fully compatible with the hydrosealt (1993) program. 
The new programs retain the same calculating options of hydrosealt plus the 
added bearing geometries, and unsteady flow and transient forced response. 

Refer to the hydroflext & hydrotran User's Manual and Tutorial for basic 
information on the analysis and instructions to run the programs. The Examples 
Handbook contains the test bearing cases along with comparisons with 
experimental data or published analytical values. 

The following major tasks were completed on 1994 (Phase II): 

a) Extension of the thermohydrodynamic analysis and development of computer 
program hydroflext to model the following bearing geometries: 

tilting-pad hydrodynamic journal bearings, 

flexure-pad cylindrical bearings (hydrostatic and hydrodynamic), and 
cylindrical pad bearings with a simple elastic matrix (ideal foil bearings). 

b) Improved thermal model including radial heat transfer through the bearing 
stator. 

c) Calculation of the unsteady bulk-flow field in fluid film bearings and the 
transient response of a point mass rotor supported on bearings. 

d) A literature review on the subject of two-phase flows and homogeneous - 
mixture flows in thin-film geometries. 
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Brief description of hydroflext program 

The motion of a cryogenic liquid on the thin film annular region of a 
fluid film bearing is described by a set of mass and momentum conservation, and 
energy transport equations for the turbulent bulk-flow velocities and pressure, 
and accompanied by thermophysical state equations for evaluation of the fluid 
material properties. Zeroth-order equations describe the fluid flow field for a 
journal static equilibrium position, while first-order (linear) equations govern 
the fluid flow for small amplitude journal center translational motions . Solution 
to the zeroth-order flow field equations provides the bearing flow rate, load 
capacity, drag torque and temperature rise. Solution to the first-order 
equations determines the rotordynamic force coefficients due to journal radial 
motions. 


The hydroflext program calculates the static load and dynamic force 
coefficients for the following bearing geometries: 

1. hydrostatic bearings with orifice compensation and rectangular recesses 

(single row or two-parallel recess row), 

2. annular pressure seals (damper seals) (cylindrical and multilobe), 

3. plain cylindrical hydrodynamic bearings (cylindrical and multilobe), 

4. fixed arc hydrodynamic bearings with arbitrary preload, 

5. tilting-pad journal bearings, 

6. flexure-pad journal bearings (hydrostatic and hydrodynamic), 

7. cylindrical pad bearings with a simple elastic matrix (ideal foil bearing). 


hydroflext includes the following thermal models: 

- adiabatic surfaces, i.e. insulated journal and bearing surfaces. 

- isothermal journal at specified temperature and insulated (adiabatic) bearing. 

- isothermal bearing at specified temperature and insulated (adiabatic) journal. 

- isothermal journal and bearing surfaces. 

- isothermal journal and radial heat flow through bearing (stator). 

- adiabatic journal and radial heat flow through bearing (stator). 

hydroflext calculates 

1) bearing flowrate or seal leakage 

2) friction torque, power dissipation and temperature rise 

3) load capacity (fluid film forces and moments) if journal eccentricity 
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is given, or 

journal equilibrium eccentricity components if the external loads are 
given. 

4) 16 complex impedance coefficients due to dynamic journal center 
displacements and journal axis rotations. The real and imaginary parts 
of the impedances correspond to the stiffness and damping coefficients 
evaluated at a specified excitation frequency. 

5) Stability indicator or whirl frequency ratio for lateral journal 
motions and equivalent stiffness at operating speed. 

6) Pressure and temperature fields on the bearing surface, and density and 
viscosity field variations, within ranges of fluid flow Reynolds 
numbers and Mach numbers. 

for 

1) isothermal flow with barotropic fluid, 

2) thermohydrodynamic adiabatic or radial heat flow and/or isothermal 
journal and bearing surfaces 

in the single phase flow regime. 

with the following fluids: 

1) liquid hydrogen, 2) liquid nitrogen, 

3) liquid oxygen, 4) liquid methane, 

5) water, 6) oil, 

7) air, 8) barotropic liquid. 

For cryogenic liquids, the fluid properties (density, viscosity and 
specific heat) are calculated with the miprops program from NBS 
Standard Reference Database 12. 

hydroflext handles the following boundary conditions at the bearing exit planes: 

1) periodic pressure asymmetry in the axial direction. 

2) local discharge end seal effects via an orifice like model to 
~~~ simulate wear-ring hydrostatic bearings or annular seals. 

3) inlet specified circumferential pre-swirl velocity distribution. 

The axial clearance functions included are of the type: 
a) uniform, b) tapered, 

c) stepped, or, d) arbitrary via spline interpolation. 

Cylindrical bearings can be specified as multi-lobe geometries, and 
bearing pads can have a specified assembly preload. 

For tilting-pads and flexure-pad bearings, pad mass moment of inertia and 
flexure rotational stiffness and damping coefficient are needed for full 
specification of the bearing geometry. 
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Brief description of hydrotran program 


The extended computer program hydrotran calculates the dynamic (transient) 
force response of a rigid rotor supported on fluid film bearings. The code 
calculates at each time step of numerical integration the unsteady bulk-flow 
field and the fluid film bearing forces due to prescribed time-varying external 
loads. The transient analysis is restricted to isothermal flows and rigid pad 
bearings, i.e., the current model does not account for transient thermal effects 
and can not handle the transient dynamics of tilting pad bearings. These 
restrictions on the code are due to the large computational times required to 
determine a transient solution with small time steps. 

The equations of motion solved are: 

M AX = FX + VX(t) (1) 

M AY = FY + WY(t) 

where M: rotor point mass [kg], 

WX,WY: external time dependent loads [N], 

2 

AX,AY(t): rotor accelerations [m/s ], and 

FX,FY(t): are the fluid film bearing reaction forces [N] . 

The external load cases considered are: 

1 periodic forcing function 

WX=Fmag cos(wt), WY=Fmag sin(wt) 

2 ramp loading along X direction 

— WX=Fmag t/T if t<T, WX=Fmag if t>=T ; WY=0 

3 step loading along X direction 

WX=Fmag for t>=0, WY=0 . 

4 rotor weight along X + sudden unbalance forcing function 

WX=Weight+ Funbal cos(«t); VY= 0 + Funbal sin(cot) 

2 

Weight=M*g , Funbal= Mrotor x UNBmag [m] x w 

5 impulse load along Y direction applied over To (secs) ^ 

WX=Weight=M*g , VY=Fimpulse exp[-.5x(t/to) ] 

The Wilson-0 method is used to perform the numerical integration of 
equations (1). A simpler program (TVOdof) uses the (linear) rotordynamic force 
coefficients at an equilibrium position and calculates an approximate transient 
response. This program executes much faster than the full numerical nonlinear 
solution and can be used efficiently for rotors supported on externally 
pressurized bearings, i.e. bearings which do not show great variation of their 
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force coefficients with the applied load or journal position. 
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ABSTRACT 

The thermal analysis of flexure pivot tilting-pad hybrid (combination 
hydrostatic - hydrodynamic) fluid film bearings for cryogenic liquid 
applications is presented. Turbulent bulk-flow, variable properties, momentum 
and energy transport equations of motion govern the flow in the thin film lands. 
Zeroth-order equations describe the fluid flow field for a journal static 
equilibrium position, while first-order linear equations govern the fluid flow 
for small amplitude journal center translational motions. Solution to the 
zeroth-order flow field equations provides the fluid film bearing flow-rate, 
load capacity, drag torque, and temperature rise. Solution to the first-order 
equations determines the linearized force coefficients due to journal radial 
motions. An elegant method to determine the tilting-pad moment coefficients from 
the force displacement coefficients, and an iterative Newton-Raphson scheme to 
balance the pad moments are detailed. Numerical predictions correlate well with 
experimental measurements for tilting-pad hydrodynamic bearings. A design of a 
flexure-pad hydrostatic bearing for a cryogenic application is shown to have a 
reducecT whirl frequency ratio without loss in load capacity or reduction in 
direct stiffness and damping coefficients. 


This report extends the previous analysis (Mid-Year Progress Report, June 
1993) of tilting- and flexure-pad bearings to include thermal effects. The 
Results section is identical to that of the original Mid-Year Report. 


For a lucid discussion on the importance and advantages of fluid film bearings 
in cryogenic turbomachinery see the Research Progress Report, Phase I, to the 
present project (San Andres, 1993). 


Refer to hydroflext User's Manual and Tutorial for information and 
instructions to run the program. Also if you are in charge of installing the 
program on your system please browse the README file enclosed on the 
hydroflext documentation. 
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INTRODUCTION 


The technology of externally pressurized bearings for cryogenic liquid 
turbopump applications has steadily advanced in the past years. Computational 
bearing design programs have been validated with a wealth of experimental data 
obtained at the TAMU Hydrostatic Bearing Test Facility (Kurtin et al., 1993, 
Childs and Hale, 1994, Franchek et al. 1994). Advanced cryogenic turbomachinery 
designs already demand an all-fluid- film-bearing technology to improve 
efficiency, reduce size and weight, and with lesser impact on the environment. 
The progress achieved brings needed verification to promote the use of external- 
ly pressurized, process-liquid bearings for use in commercial high-performance 
rotating machinery. This type of machinery eliminates expensive oil-based 
lubrication systems as well as rolling element bearings of limited success in 
high-speed applications. Magnetically driven canned and sealess pumps of proven 
superior performance and reliability are examples of these new advances. In 
other high performance applications, hydrostatic bearings and damping bearing 
seals are the primary candidates to provide this quantum step in machinery 
performance. 

Despite the attactive features provided by hybrid (combination hydro - 
static and hydrodynamic) journal bearings (HJBs), fluid film bearing stability 
considerations and thermal phenomena accompanied by fluid phase change are 
issues of primary concern for high speed operation with large pressure 
differentials (San Andres, 1993). Bearing dynamic stability is related to 
hydrodynamic and liquid compressibility effects limiting severely the maximum 
operating speed of a rotating structure to a value close to the system first 
critical speed divided by the whirl frequency ratio (VFR) . The WFR is a 
stability indicator which, for conventional incompressible liquid hybrid 
bearings, is typically equal to 0.50. Innovative fixed geometry hybrid bearings 
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with improved stability have been designed and tested (Franchek, 1992). Some of 
these bearing configurations include rough-surface bearings, anti-swirl liquid 
injection HJBs, and bearings with structural stiffness asymmetry. 

The flexure-pad hybrid bearing concept (Zeidan, 1992), see Figure 1, 
constitutes a promising alternative for use in high speed applications and 
without stability restrictions. The major advantage of this bearing over 
conventional tilting-pad technology derives from its structured design and 
single-piece EDM manufacturing. The thin webs provide enough radial stiffness to 
support aplied loads while allowing for rotational pad motions. Tilting-pad 
bearings are mechanically complex due to their many parts, and on assembly, 
tolerances stack-up requiring a high degree of manufacturing accuracy to yield 
acceptable clearances. Furthermore, these bearings are supported on ball and 
socket type joints which are prone to rapid wear if the fluid media has no 
lubricity. This is precisely the case with liquid cryogens as well as with most 
process liquid applications. The flexure-pad structural arrangement offers a 
unique advantage to ensue combined hydrostatic and hydrodynamic actions with 
supply ports directly machined on the thin webs for delivery of a pressurized 
lubricant. The flexure-pad bearing technology applies advanced automated 
machining processes to manufacture bearings of any specified tolerance for shape 
and film clearance as well as structural properties (Zeidan, 1992). In addition, 
the dynamic variation of external pressure and flow through the bearing orifice 
ports can lead to the concept of a SMART BEARING, i.e. a self-adapting, low 
friction, mechanical support able to modify the entire rotor-bearing system 
frequency response by controlling its stiffness and suppressing critical speeds. 
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Literature Review 


Tilting-pad hydrodynamic journal bearings are used frequently as support 
elements in high speed rotating machinery requiring maximum stability. The 
ability of the bearing pads to rotate (about an individual pad pivot) determines 
a reduction or complete elimination of the undesirable cross-coupled stiffnesses 
which produce bearing instability in fixed geometry fluid film bearings. Flack 
et al. (1984), and Fillon et al. (1987, 1991) offer lucid reviews of the 
extensive theoretical literature relevant to tilting-pad bearings. Shapiro et 
al. (1977) and Nicholas et al.(1979, 1994) discuss practical design issues in 
tilting-pad bearings and refer extensively to the benefits these bearings have 
brought to rotating machinery. 

The theoretical models used to describe the performance of tilting-pad 
journal bearings have steadily increased in complexity. The number of 
independent parameters required to fully represent the flow are virtually 
countless. Three-dimensional thermohydrodynamic (THD) and structural effects as 
well as the lubricant feed and mixing mechanism are thought to be necessary for 
accurate prediction of peak pressure, minimum film thickness and temperature 
rise. Theoretical models currently account for variable properties, heat 
transfer through the bearing pads, pad thermoelastic deformations, and pivot 
wear and flexibility. The analyses of Ettles et al. (1980, 1992), Knight and 
Barrett (1988), and Taniguchi et al. (1990) are representative of these advanced 
developments. Pinkus (1990) offers a scholarly account of the major developments 
in the area of thermal effects in fluid film lubrication. 

Tilting-pad bearings have many more degrees of freedom than fixed 
geometry bearings, and consequently, their analysis is complicated and subject 
to a number of assumptions regarding the dynamic motion of journal and 
pads. Lund (1964) introduced rotordynamic force coefficients for tilting-pad 
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bearings using an assembly method from a single generic fixed pad. Barrett et 
al. (1988) and White and Chan (1992) discuss the effect of excitation frequency 
on the bearing reduced force coefficients. Rouch (1983) includes pivot 
flexibility, and Lund and Bo-Pedersen (1987) and Kim and Palazzolo (1993) 
account for pad flexibility on the THD evaluation of rotordynamic coefficients. 
These investigations show that thermoelastic deformations affect significantly 
the bearing damping force coefficients. Nicholas and Barrett (1986) also 
demonstrate, both theoretically and experimentally, that support flexibility has 
a profound effect on the critical speeds and dynamic performance of rotors 
supported on tilting-pad bearings. 

Orcutt (1967) and Pinkus (1984/85) present analyses for prediction of 
tilting-pad bearing performance for turbulent flow and isothermal conditions. 
Taniguchi et al. (1990) provide tilting-pad bearing pressure and temperature 
experimental data and detailed correlations to numerical predictions from a 
complex three-dimensional, turbulent flow, thermohydrodynamic analysis including 
heat conduction through the bearing pads. Bouard et al. (1994) compare different 
turbulent flow models to the test data of Taniguchi et al. with no major 
diffrences among the models used. Recently, Ha et al. (1994) include a pad 
leading edge-ram pressure effect to improve the prediction of pressure and film 
thickness on large scale, turbulent flow tilting-pad bearings. 

Relevant experimental investigations in tilting-pad bearings are given 
by Someya (1988), Brockwell and Leinbub (1989), and Fillon et al. (1991) for 
laminar flow conditions, and by Orcutt (1967), Taniguchi et al. (1990), and Ha 
et al. (1994) for both laminar and turbulent flows. The detailed experimental 
works have concentrated on the static performance of tilting pad bearings under 
various conditions of load and shaft speed with attention to the film and pads 
surface temperature and pressure evolution. In general, theoretical predictions 
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agree well with measurements if all bearing parameters are known accurately. 
Actual operating clearance and preload, and the inlet pad temperature are 
perhaps the most important factors affecting tilting-pad bearing performance. 
Isothermal models suffice to predict well bearing performance for small loads 
(or large Sommerfeld numbers), while apparently a complete elasto-thermohydro- 
dynamic model becomes mandatory for large loads (small Sommerfeld numbers) 
accompanied by very small minimum film thickness and large temperature rises. 

Measurements of tilting-pad journal bearing force coefficients have been 
relatively scarce in the technical literature until recently, see for example 
the relevant works of Someya (1988), Brockwell et al.(1990), Parkins and Horner 
(1992), and Arumugan et al. (1994). Test stiffness coefficients agree reasonably 
well with theoretical predictions. Measured damping coefficients are typically 
lower than calculated values. Pad thermoelastic deformations, and pivot and 
shaft flexibilities appear to affect greatly the measured results. Flack and 
Zuck (1988) report experimental rotordynamic responses of two flexible rotors 
supported on tilting-pad bearings. These fundamental tests demonstrate unstable 
rotor -be aring behavior as a whip condition with whirl frequency equal to the 
first critical speed. The instabilities are thought to be bearing induced, but 
the physical mechanism or bearing parameter(s) producing this highly unusual 
behavior (in tilting-pad bearings) are yet to be identified. 

Zeidan (1992), and Zeidan and Paquette (1994) introduce flexure-pivot, 
tilting-pad bearings and discuss their advantages over conventional tilting-pad 
bearings. Armentrout and Paquette (1993) and Chen (1994) present analyses for 
evaluation of rotordynamic force coefficients in isothermal flow, flexure-pad 
bearings. Experimental bearing force coefficients are yet to be reported in the 
archival literature. Zeidan (1992) and Chen (1994) describe applications in 
which flexure-pad bearings provided stability to otherwise unstable rotating 
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machinery and with a performance superior to that of conventional tilting-pad 
bearings. De Choudhury et al. 1992) demonstrate experimentally that a flexure- 
pad bearing brings a lower temperature rise and frictional power loss when 
compared to a similar size five shoe tilting pad journal bearing. This fact 
allowed the flexure-pad bearing to operate at a lower oil flow rate while 
providing acceptable oil throwoff temperatures. 

The analytical work on fluid film bearings at Texas A&M University has 
lead both the experimental verification and the practical application of hybrid 
bearing technology. Many of the relevant and some times surprising analytical 
results have been later verified on the test stand. Most notably the reduced 
whirl frequency ratio (=0.50) has been of primary concern (San Andres, 1990). 

The measures to fix or avoid the instability condition which limits severely the 
application of hybrid bearings to cryogenic turbopumps have been discussed in 
detail elsewhere (San Andres, 1993). Tilting-pad journal bearings with 
hydrostatic pressurization were originally proposed by Greenhill (1991) as ideal 
bearing candidates with no hydrodynamic stability restriction. However, the 
mechanical complexity of conventional tilting pad bearings prevented their 
consideration as a suitable replacement for fixed geometry hybrid bearings. The 
appearance of flexure-pivot tilting-pad bearings on the commercial bearing 
market has relaxed the concerns for their potential use in cryogenic liquid 
turbopumps. The present technical development then presents the relevant 
analysis and computational predictions which demonstrate the feasibility of this 
type of bearing configuration for aerospace applications. 
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ANALYSIS 


Consider the flow of a variable properties liquid in the thin film annular 
region between an inner rotating journal and a hydrostatic/hydrodynamic tilting 
pad bearing. The bearing pads provide rotational stiffness and viscous damping 
(see Figure 1) but are regarded as rigid in the radial direction. Each pad may 
contain one or more recesses and orifice restrictors for radial liquid supply to 
the bearing. A liquid at high pressure (P g ) is supplied radially through 
orifices and flows into the bearing recesses. For hydrostatic operation, the 
recess pressure (P ) is uniform over the entire recess volume, while at the 
recess edges, an inertial pressure drop occurs due to the sudden transition from 
the recess of depth (H^) into the film land region of small thickness (H). For 
hybrid operation (combined hydrodynamic - hydrostatic), the pressure within the 
recess downstream of the supply line rises due to viscous shear effects. The 
viscous pressure rise occurs for the small recess depths typical of cryogenic 
HJBs and closely parallels the pressure development in step bearings 
(Constantinescu et al., 1975, 1987). The liquid then flows into the film lands 
and a~pressure drop to the discharge (sump) value follows. On the thin film 
region, hydrostatic and hydrodynamic effects are important, and for the high 
speed - high pressure conditions typical of cryogenic liquid environments, the 
fluid flow in the bearing is fully turbulent and dominated by fluid inertia 
effects. 

The analysis studies the fully developed turbulent bulk-flow of a cryogenic 
liquid whose material properties depend on its local thermophysical state of 
pressure and temperature. Cryogenic liquids have low viscosities, and 
thermal (energy transport) effects due to friction heating and kinetic energy 
variations are expected to be of minor importance for pure hydrostatic 
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performance. This assertion is not fully justified for especial operating 
conditions (see for example San Andres et al., 1993, and Yang et al., 1993). 
Coordinate system and film thickness 

A local coordinate system is placed on the unwrapped plane of the bearing 
with the {x,z} axes pointing in the circumferential and axial directions, 
respectively. This cartesian coordinate system is a logical consequence of the 
smallness of the film thickness relative to the bearing diameter and length, and 
constitutes the fundamental assumption of lubrication theory. 

A characteristic tilting pad of angular extent 8^ and axial length equal to 
L=V l r is shown in Figure 2. The leading edge of the k-th pad is defined by the 
coordinate 0^^, and 9^ denotes the position of the pad pivot or rotational 
flexural node . At operating conditions, the journal position relative to the 
bearing housing is described with reference to the inertial axes {X,Y,Z} by the 
journal center displacements (e x ,e Y ). Simple geometrical relationships determine 
the film thickness in the flow region of the k-th bearing pad to be given by: 

H (Z,97t)= Cp(Z) + cos(0) + e Y sin(0) - r^cos^- 6 P ) - R* S^sin( 0- 6p) (1) 

where 0=x/R and 6 is the pad rotational angle; is the pad machined radial 

clearance, in general a function of the axial coordinate (Z); r =C -C , 

ppm’ 

corresponds to the bearing preload with C m known as the bearing assembled 
clearance. The journal center displacement components (e x ,e Y ) and pad rotational 
angle (8 ) are functions only of time. Note that no journal misalignment has 
been considered on the description of the film thickness. This ommission is due 
only for simplicity in the presentation of the analysis. 
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The equations of motion 


The fluid flow on the film lands of a bearing pad is considered as fully 
developed with a turbulent character due to the large axial pressure drop, high 
rotor surface speed and the low viscosity typical of cryogenic liquids. The 
equations of motion for the turbulent bulk-flow on the thin film land of the 
k-th bearing pad are given by (San Andres, 1992, 1993, Yang, 1992): 


equation of continuity: 


a a 

— (pH) + — 


3t 


(PHU ) k + (pHU fl ) k = 0 

3z Z 3x 9 


( 2 ) 


axial momentum equation : 


^ ( K „ l\ iKl,H1 V k | ( 8( ‘ ,m W + 3 <° HU 8 U Z > f 

3z H k l Z Z / 3t 1 3z 3x j 


(3) 


circumferential momentum equation: 


3P k y 


k M | m k 3(pHU 0 ) k ( 3(pHU z U 0 ) 3(pHU 0 U 0 ) ) k 

-H = — , < K . > + + + > ( 4 ) 

3 x H k l 0 9 J 2 J 3 t l 3 z 3 x J 


energy tranport equation: 


l3t 


3 ( ' 

- pHT 

V / 


3x . 

l 


\ 


pHU.T 


(3P 


2 2 1 
U +U + =-0-R*U 
x y 2 x 


+ KjQR 


3P ) k 

H 3P 

— + 

Q * R 

3x. J 

2 3x 

1 


f N 

4 

1 k 

j-QR-U 
4 x 

1 , (5) 


i=x,y 


k=l , . . . ,N 


pad 


on the flow region {-L^<Z<L^; ; and where, K z =K 0 =(kj+kg)/2 are 
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the wall shear stress difference coefficients as local functions of the 
turbulent friction factors, bearing and journal surface conditions, and the flow 
Reynolds numbers relative to the rotating (Rj) and stationary (R g ) surfaces, i.e 

1973). The cryogenic liquid properties are calculated 
from the Benedict-Web-Rubin equation of state as given in the standard data base 
of McCar ty(1986) . In the energy equation, 

V h B (T - T B> + h J < T - T J> (6) 

denotes the heat flow from the fluid film to the bearing and journal surfaces at 
temperatures (Tg) and (Tj), respectively, and with (h^hj) as the bulk-flow 
convection heat transfer coefficients. A description of these parameters in 
turbulent bulk-flows is detailed by Yang (1992). 

Hydrostatic Bearing, Recess Floy and Pressure equations 

A mass conservation equation at each bearing recess of area (1*R*0 ) k and 
depth is defined by the global balance between the mass flow through the 
orifice restrictor (Q ro )> the mass flow into the film lands and the time rate of 
change of liquid mass within the recess volume V^_. This equation is given as; 
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2 

where A Q =C d nd /4 is the effective orifice area, and e p =(l/p)3p/3P and 

1/ p) 3 p/ 3T represent the fluid compressibility and volumetric expansion 
coefficients, respectively. represents the closure of the recess volume with 
the film lands and with normal n along the boundary line. Note that the orifice 
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flow equation is valid only for small changes of the liquid density. A more 
accurate relationship should be used for large variations in fluid material 
properties across the orifice. A fundamental relationship is yet unknown for 
compressible cryogenic liquids, although a general thermophysical model is 
available in the literature (Hall et al., 1986). In equation (7) above, the 
orifice discharge coefficient is of extreme importance since experimental 
measurements at Texas A&M University and elsewhere have shown that numerical 
predictions require accurate values to provide meaningful results. Needless 
to say that this orifice coefficient depends not only on the flow structure and 
Reynolds numbers, but also on the fluid used and the closeness of the orifice to 
the journal surface. Refer to Kurtin et al. (1993), Franchek (1992) and Mosher 
(1992) for experimentally based discharge coefficients in high-pressure, 
high-speed water HJBs, and to Scharrer (1990) for some measurements on liquid 
freon and liquid nitrogen HJBs. 

The fluid edge pressure at the entrance to the film lands is given by the 
superposition of viscous shear effects on the recess extent and an entrance drop 
due to~*Tluid inertia. Figure 3 shows the assumed pressure distribution within 
the recess volume and details the relevant nomenclature. On the circumferential 
direction, the pressure rise (P~) downstream of the recess orifice is given by 
(Constantinescu et al., 1975, 1987, San Andres, 1992): 
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where, M is the circumferential flow local Mach number at the orifice discharge 
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and defined as the ratio between the azimuthal velocity and the liquid sound 
speed (l.//[f3 p p r ]). 

The entrance pressures (P^) to the film lands in the circumferential and 
axial directions are given by: 

P e - % - “e 

■ p r - f (»:/p;) 2 n 2 } u 2 

for r=l,2..., Nrecess on k-th pad 

The analysis generalizes equations (8) for uneven empirical entrance loss 
factors £ in the axial direction (Z) and also circumferentially upstream(u) and 
downstream(d) of the recess. The Bernoulli like pressure drop occurs only if 
fluid flows from the recess towards the film lands. On the contrary, if fluid 
enters from the film lands into the bearing recess, then the edge pressure is 
regarded as equal to the recess pressure. 

A simplified form of the energy transport equation at a bearing recess is 
esily derived from equation (5). A detailed exposition of the energy transport 
process in a bearing recess is given by San Andres (1993). 

Boundary Conditions 

At low rotational speeds, the pressure at the side (right and left) 
boundaries of the k-th pad is essentially constant and equal to a specified 
value of ambient or sump pressure, i.e., 

P k (L R ,9) = P k ; P k (-L L ,9) = P k (9. a) 

and at the pad leading and trailing edges, 

P k (Z,^ k ) . P k <Z,<W - f, <P Ei P L ) on -L l < Z < Lg (».b) 
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where (f s ) denotes a linear pressure variation along the axial side of the 
bearing pad leading and trailing edges. At high journal surface speeds (Q-R) 
significant momentum changes occur at the pad leading edge (0^). Inmediately 
upstream of the pad, the fluid entering the film lands can develop a dynamic 
head equal to some fraction of a reference dynamic pressure based on the bearing 

surface speed (Burton and Carper, 1967, Smalley et al., 1974, Mori et al., 1991, 
Ettles et al., 1968), i.e. 


4plC(z ’ 9 ik>- - y«>- 


ra-R) 


l 2 


(9.c) 


The coefficient (K p ) is an empirical ram pressure factor. Burton and Carper 
(1967) suggest a value of K p =0.64 for high speed flows with large turbulence 
levels. The appearance of this ram pressure effect due to journal rotation at 
the leading edge of a bearing pad is of fundamental importance on the analysis 
of high speed pad bearings. The model introduced is very simple. A more accurate 
knowledge of the entrance entrance condition requires the analysis of the 

complex flow field on the axial groove and feeding port separating the bearing 
pads (Ha et al. , 1994). 


A further assumption about the external flow is necessary to fully specify 
the flow variables. The flow at the inlet(leading edge) to the film land in the 
pad bearing is assumed to be parallel to the direction of surface motion. This 
simplification implies that on the line, 


- L L^%’ U z< Z -«lk> - 0 . k.l,..,N pad (10) 

In the thermal analysis, the inlet bulk-flow temperature to a bearing pad 
is obtained using mass flow and energy balances of the upstream pad flow and 
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temperature and incoming fresh fluid at a supply temperature (T s ). The pad inlet 
bulk-flow temperature depends on the specification of a thermal mixing parameter 
as given by Mitsui et al. (1983). 

The boundary conditions noted refer only to the fluid film flow within the 
bearing pad. At an equilibrium journal position due to an imposed external load, 
each individual tilting- or flexure-pad must be balanced, i.e. the fluid induced 
moments and elastic restoring moments must add to zero. This item will be 
discussed later in the analysis. 

Perturbation Analysis 

Consider the motion of the journal as the superposition of small ampitude 
periodic motions around an equilibrium static position. That is, we let the 
journal center displacements and pad rotation be given as 


e X<*> = e Xo + Ae X el<0t ’ V” " e Yo + 
S k (t) = S k + A6 k e iwt , uih 


( 11 ) 


where denotes the frequency of the whirl motion. The magnitudes of the dynamic 

perturbations in journal displacements and pad rotations are very small, i.e., 
k, 

)|/C* <<< 1. Then, the film thickness can be thought as the 
superposition of steady-s tate(h Q ) and dynamic (h^) components given by the real 
part of the following dimensionless expression: 
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or 


with 
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where C* is a characteristic clearance for the bearing. The functions defined in 

equations (14-15) facilitate the understanding of the perturbed flow field 

equations and the resulting rotordynamic force coefficients given latter. 

For each bearing pad the flow field variables (U ,U~,P,T), as well as the 

z y 

fluid properties (p,u) and the shear parameters (K~ f K ) are also formulated as 

y z 

the superposition of zeroth-order and first-order complex number fields 
describing the equilibrium condition for steady-state flow, and the perturbed 
condition for small amplitude dynamic journal motions, respectively. In general, 
these fields are expressed as: 


Y k = e 1Wt |Ae x i' x +Ae Y '{' Y +(R/C^)ASY s j= + e 1Wt A a (16) 

a=X,Y,8 k , k=l,2..N pad 

Substitution of (12) and (16) into the equations of motion (2-7) determines the 
zeroth-order dimensionless equations : 
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k=l,2..,N . 

’ pad 


Refer to the Nomenclature for definition of all dimensionless variables in 
equations (17-20). On the recess volumes of a hydrostatic bearing , the mass flow 
through the orifices must equal the mass flow into the bearing lands, 
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r=l,2.. Nrecess on k-th pad 

The circumferential pressure rise at the downstream side edge of a recess is 
related to the recess pressure and fluid velocity by: 
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where ‘f'j.-K^'S^/^'h^) is a dimensionless recess shear factor. The recess edge - 
entrance film pressures in the circumferential and axial directions are 
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given by 


ik. 


+ - e. 2 

^eo ^eo ^ U B° 


nk 


+ p e , 2 

^eo ~ ^ro ^ ^ z Uzo 


where, 


*e = (1+ V' Re P 


i-(Pg/p;)n 2 


1.95 


1 +- 


Re, 


0.43 


and 


+z = ( 1 + V' Re p 


1_(p e /p e )n2 


(22. b) 


(22. c) 


(22. d) 


(22. e) 


At the bearing pad discharge (axial) sides and the trailing edge, the 
pressure takes specified values, see equations (9). The ram pressure jump at the 
leading edge of a bearing pad is given in dimensionless form as: 
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Across the line, -L^/R< z <L^/R, (pad leading edge) the inlet flow is regarded 
as parallel to the bearing pad, i.e., 

u zo (z ’ e ik ) = 0 (23 - b) 

First-order dimensionless equations 

The first-order flow equations for perturbations in the film thickness due 
to journal displacements (Ae^,Aey) and pad rotation (A5^) can be expressed in 
the unique form: 
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for perturbations in the a=X,Y directions and pad rotation oc=S k , k=l, 2 ..,N 

’ ’ pad 

In the equations above (h a ) correspond to the perturbation fields in film 
thickness as defined by equations ( 15 ). The coefficients y^j , s arise from 
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perturbation of the shear stress factors k z , k 0 and k Jf and their explicit 
formulae have been given before (San Andres, 1993, Yang, 1992). The perturbed 
^* rs * _or ^ er fluid peoperties are obtained from relationships like: 

? a = O//3p)| 0 p a + 0?/3f)| o f a ; /=p,u,c p (28) 

On each recess of a kth-pad hydrostatic bearing, the first-order 
continuity equation takes the linearized form: 


r ,2 - 

* p ro 


ra 


-ik 


L 2q 


1 P Po (1_P ro ) [ +i ' p ro X 2 


ro 


■ \ 

= 1 • p A,,' 

ro 11 


h d9 + 
ra 


Ir 


f 


P h u +p h u +p u h Ln dT + 
o a o o o a a o o [ r 


(i-p ) 
l ro / 


S p 2_ 
* ro 


2q 


- 1X lP; 


ro 


ro 


»ro T 


ra 


r=l,2. . .Nrecess on k-th pad 


(29) 


The first order recess edge (circumferential and axial) pressures into the film 


lands are given as 
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for a=X,Y,5 k 

No perturbations on the pad leading edge ram pressure are considered. Thus, 
the first-order pressure field is null everywhere on the side boundaries of the 
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bearing kth-pad; i.e. p k =0. 

a 

Note that the first order equations, (24) to (30) can be written in the 
general form: 
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where L are linear differential operators with coefficients depending only on 
the zeroth order flow field variables (u^, u 2Q , P o ,h Q ) k . Equation (31. c) 
describing the perturbed flow field due to a pad rotation (AS k ) is just a linear 
combination of the first two equations. Hence, for homogeneous boundary 
conditions on the perturbed fields, the following statement is correct: 


= sln(8 P ) ' i x - cos( 0 d> ' l r where V u es’ u z5’ p 


(32) 


This equation shows that the first-order flow field due to a pad rotation can be 
found as the superposition of the perturbed flow fiels due to journal 
displacements in the {X,Y} directions. Thus, the computational time for the 
evaluation of force coefficients in tilting (or flexure) pad bearings is 
considerably reduced. 

^ Fluid Film Forces, and Dynamic Force and Moment Coefficients 

The fluid film forces and moment on each bearing pad are calculated by 
integration of the pressure field, i.e., 
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where P is a characteristic pressure. The pad moment due to fluid film effects 
sa 

can be rewritten as: 


M k = R | sin(e£) F k - cos(fe p ) F k j 


(34. b) 


Substitution of the pressure field, equation (16), 
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into equations (33) and (34) determines the k-th pad zeroth- and first-order 
forces and moment due to the fluid film pressure as: 
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and 
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are the zeroth order fluid film forces and induced moment for the kth-pad, 
respectively* The are pad impedance coefficients whose real and imaginary 
parts correspond to pad fluid film stiffness and damping coefficients, 
respectively* Analytical expressions for the kth-pad impedance coefficients are 
given as: 


force impedances due to journal displacements Ae v , Ae v : 
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force impedances due to pad rotation AS * 
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pad moment imp edances due to journal displacements Ae v , AeY: 
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A total of 9 impedance coefficients are determined for each bearing pad. A clear 
distinction of fluid inertia force or added mass coefficients is not appropriate 
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for tilting-pad bearings since the stiffness and damping coefficients are 
frequency dependent. A major simplification in the evaluation of the impedance 
coefficients is possible by using equation (32) for the dynamic pressure field 
due to a pad rotation, i.e 

p 6 = sin <®£> p x ~ Py> 


and allows the impedance coefficients to be expressed as unique functions of the 
force impedances Z^, a, 0=X,Y as follows: 

k k k k k 
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Equations (39) show only four impedance coefficients to be independent, with the 
other five uniquely determined by the relations defined above. These equations 
simplify greatly the analysis of tilting-pad bearings since only the dynamic 
pressure fields for journal center perturbations in the (X,Y) directions are 
(strictly) required. 
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Bearing Equilibrium Equation and Pad Equation of Rotational Motion 
In a multiple pad bearing, the sum of the individual pad forces must balance the 
externally applied load V. This load has components (V^,Wy) on the {X,Y} 
directions, respectively. Thus, the static equilibrium force equation (in the 
absence of rotor-inertial effects) is given as: 
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(40) 


The equation of rotational motion for each bearing pad is given as: 

*k ^ * Cr k* k * Kr k sk ’ R I sinOp) f£ - cos(6p) P R ) (41) 

k.l , . . . N pad 

where 1^, Cr^, Kr^ correspond to the pad inertia, and pad rotational viscous 
damping and structural stiffness coefficients, respectively. (Note that in a 
conventional tilting pad bearing, Kr,=0). 

1C 

A static load W q impossed on the fluid film bearing determines a bearing 

equilibrium position at a static journal eccentricity (e v , e v ) and a set of 

xo Yo 

stationary pad rotations {S Q } ,1 c =l,..Np a( j. Hence, the bearing force and pad 
moments equilibrium equations are givens as: 
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0 ■ _Kr k^ * ■ - Kc k 5 o * R ( sin <^> Pjo" cos < e ^) f Yo) < 43 > 

k=l , 2. . . , N pad 

In the practical analysis of tilting- and flexure-pad bearings, reduced bearing 
stiffness and damping coefficients are determined at the static equilibrium 
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position and for a particular frequency of excitation, typically synchronous 
(«=Q) or subsynchronous (w<Q). For small amplitude journal motions (Ae Y , 

Aey) e , and pad rotations k=l, N ^ with frequency c*> about the 

equilibrium position (e Xo , the individual pad forces and the pad moment 

balance equation can be written as (see also equations (35)): 
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where + i » are the pad impedance coefficients at frequency «. 

Equation (45) indicates that the kth-pad rotational motion is a function of the 
journal dynamic displacements (Ae^jAey)- 


Let 


ZZ SS * (K 55 tKr k-“ 2l k> + 1 “ < C S5* Cr k> 


(46) 


and the relation between the k-th pad dynamic rotations and journal dynamic 
displacements is given from equation (45) as: 
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Substitution of (47) into the force equation (44) renders the pad reduced 
impedance coefficients as: 
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a,P=X,Y (49) 


and > <W ^RaP^ a, P=X,Y are re duced force-displacement stiffness and 
damping coefficients for the kth-pad. 

The reduced stiffness and damping coefficients for the entire bearing 
assembly are then given by the summation of all the pad reduced force 
coefficients, i.e. 


'Stop " & K R«P 5 Slap = |!l d C R«P «»3=X,Y (50) 

A method to determine the bearing static equilibrium position 

A static load with components (V Xo ,V Y q ) determines the journal center 
euilibrium coordinates (e Xo .e yo ) and pad rotation angles {5 k }, k-1,.. N pa(J ,such 
that equations (42) and (43) are satisfied. In practice, an iterative method is 
needed to determine the equilibrium journal coordinates and pad rotation angles 
satisfying the applied load and pad moment constraints. The current analysis 
implements a Newton-Raphson algorithm to reach convergence to the actual 
solution. Consider at iteration J, the journal center to have an ecentricity 
defined by its components ( e x > e y) on the { X , Y] directions, and a set of pad 

k J 

rotational angles {5 } , k=l,..,N ^ such that the pad moment equation (43) is 

satisfied, i.e. the pad moments are balanced. For these conditions, the bearing 
has a set of reduced static stiffness coefficients (calculated at o>=0) denoted 
aS {K Rap )J ’ a ’P= X ’ Y - 

Improved journal center coordinates (e x ,e y ) J+1 converging to the actual 
values satisfying equation (42) are determined as: 


J+l J a J+l J 
X X + ^X 5 e Y = e y + ^Y 


(51) 
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where the small 
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Improved pad rotations are found from: 
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N °te that at each step in the iterative process the pad moment equation is 
satisfied, i.e. , 
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The procedure continues with updated journal eccentricities and pad angles until 
the error in the load equation is within some predetermined range, i.e. 

l W Xo +F X I - e tf 0 ’ l U Yo +F Y J | - e W o where e=0-01 typically. 

The load-bearing equilibrium position scheme described here has proven to 
be robust because at every iterative step each pad moment equation is satisfied, 
and the improved values of journal eccentricity and pad rotation angles are 
predicted from the stifness coefficients evaluated at zero whirl frequency. 

It is important to note that most algorithms for tilting-pad bearing analysis 
sometimes have difficulty in determining the journal equilibrium position when 
one or more pads become unloaded (*). This is most frequent for bearings with 


J ? f A pad j S deflned as unloaded (fully cavitated) when no fluid pressures 
(also forces and moments) are generated within the pad for the values of journal 

^ ^ C i ty i? nd £ ad r ° tatlon an « le specified. At this condition, it should be 
clear that all pad stiffness coefficients are zero, in particular K 

oo 
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small, null or even negative preload. The traditional approach artificially 
increases the pad rotational angle until this becomes "barely" loaded. This 
procedure only works for tilting-pad bearings, i.e. pads without rotational 
structural stiffness. On the other hand, for flexure-pad bearings, if the pad 
becomes unloaded (null fluid film moment), then the only admissible equilibrium 
solution is S =0. 

Another more difficult problem arises when a pad goes from unloaded to 
loaded in the current iterative step. Here, unless the predicted (new) pad 
rotation angle is sufficiently large, the iterative scheme will bring the pad to 
a state of unloading when the actual solution indicates the opposite. The 
procedure then introduces spurious oscillations, and frequently takes large 
computing times to reach an "erroneous" solution. The difficulty described is 
overcome easily by recognizing that a valid equilibrium (stable) pad rotational 
angle requires a positive pad moment stiffness, i.e. 

k V K S5 > 0 < 55 > 

Numerical method of solution for the flow equations 

The well known control-volume method of Patankar (1980) is used to solve 
the differential equations of motion. Staggered grids containing control volumes 
for the primitive flow variables: circumferential and axial velocity, and 
pressure and temperature are generated. Algebraic difference equations are 
derived on each control volume for the conservation of circumferential and axial 
momentum, and transport of energy. A discrete form of the conservation of mass 
equation allows determination of the pressure field. The SIMPLEC procedure of 
Van Doormaal and Raithby (1984) is adopted for solution of the non-linear 


37 



difference equations along with a scheme for line solution developed by Launder 
and Leschziner (1978). Further details on the method applied for hybrid 
(hydrostatic and hydrodynamic) bearing configurations are given by Yang (1992). 

Figures 4 shows the thermal bulk-flow models implemented on the 
computational program. These models include adiabatic (insulated) journal 
surface or specified (isothermal) journal temperature, as well as isothermal 
bearing surface or bearing shell with radial heat flow and including an 
adiabatic surface. Complete details on thermal boundary conditions and numerical 
treatment are given by San Andres (1993). 

The numerical scheme devised to find the equilibrium solution in tilting 
pads and flexure pad bearings is both elegant and robust. Non-reduced force 
coefficients for each bearing pad are found from perturbed analytical difference 
equations for journal displacements in two orthogonal directions. The pad moment 
coefficients and force coefficients due to pad angular motions are fully defined 
as linear combinations of the force coefficients due to journal dynamic 
displacements, see equations (39). The method then avoids numerical evaluation 
of monism force coefficients due to pad rotations and thus, it is fast and 
accurate. Recognition that the overall pad moment-rotation stiffness coefficient 
has to be positive for a true pad equilibrium solution brings robustness to the 

procedure and avoids the possibility of oscillatory solutions, in particular for 
unloaded pads. 


RESULTS and DISCUSSION 

The validaty and usefulness of the analysis and computational program 
developed can only be assesed by extensive comparisons with relevant experiment- 
al data. In the absence of test results, comparisons with predictions from 
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similar analyses given in the literature are then justified. In any case, the 
calculations from the computational model implemented can only be as good as the 
assumptions used in the analysis. It is noted here that the present theoretical 
development has as its major objective the evaluation of the static and dynamic 
force performance characteristics of cryogenic liquid fluid film bearings. This 
particular application is known to require large journal surface speeds or 
external pressurization to provide the required load capacity since cryogenic 
liquids are essentially inviscid. At this time there are no experimental results 
for til ting-pad or flexure-pad hybrid bearings handling cryogens. Measurements 
of forces and force coefficients for fixed-arc (i.e. rigid surface) hybrid 
bearings with cryogens and surrogate fluids have appeared recently. The 
experimental work of Childs and Hale (1994) provides a wealth of test results 
and pointed out new directives for research and uses in the aerospace industry. 
Previous publications from the author, San Andres et al. (1990-1993), present 
extensive comparisons with measurements obtained from the hydrostatic test 
facility of Childs, as well as from other relevant investigations (Butner and 
MurphyTl986, Adams et al., 1992, Chaomleffel et al., 1986). 

The experimental work of Taniguchi et al. (1990) is most relevant to the 
present development since the tilting-pad bearing investigated operates both in 
the laminar and turbulent flow regimes. Furthermore, the measured pressure shows 
a "ram" effect due to fluid inertia at the leading edge of the loaded bearing 
pads. The THD analysis of Taniguchi et al. accounts for the effects of viscosity 
variation across the film thickness, includes a turbulent flow model based on 
the Ng and Pan Theory (1969), and a three dimensional energy tranport equation 
for the fluid film with a differential heat conduction equation for the bearing 
metal. The model imposses a realistic boundary condition at the shaft interface 
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and allows for a thermal mixing flow model at the pad leading edges. The 

analysis does not include thermoelastic pad deformations. The importance of this 

effect on tilting pad bearing performance has been clearly ellucidated by Fillon 
et al. (1991). 

Table 1 presents the geometry and operating characteristics of Taniguchi 
1990) four shoe tilting pad bearing. The large bearing tested (19 inches 
diameter) is for use in a steam turbine and was loaded to 300 kN (67,500 lbs). 
Numerical calculations using the present model have been performed using the 
data shown in Table 1 with an adiabatic bearing and shaft surfaces thermal 
model, fluid inertia at the film lands and a pad inertia entrance factor 
K p -0.25. A pad heat carry-over mixing coefficient equal to 1.0 was assumed in 
all the calculations. Figure 5 shows the equilibrium journal eccentricity as the 
load increases from 60kN to 300kA for a shaft speed of 3,000 rpm. The predicted 
results agree very well with the measurements and also with the theoretical 
results from Taniguchi et al. (Figure 8 of referred paper). Figure 6 shows the 
calculated and experimental minimum film thickness for the loaded pads (2 and 3) 
for th-e-same range of applied external loads. The agreement between the present 

numerical results and test data follows the same trend as the theoretical values 
provided by Taniguchi et al. 

At a fixed load of 180kN and shaft speed of 3,000 rpm, Taniguchi et al. 
present detailed measurements of the bearing centerline pressure, film thickness 
and metal surface temperature. Measured frictional drag loss is equal to 200klt 
and the present numerical calculations predicts a similar value. The test 
results correspond to a (slightly, turbulent flow case with a circumferential 
flow Reynolds number (p2RC/u> equal to 1,430. Figures 7 through 9 show 
predictions of pressure, film thickness and bulk-flow film temperature from the 
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present model for two cases. The first case labeled as conventional shows 
results obtained with a model not including fluid inertia effects at the film 
lands or at the pads leading edges. This case corresponds to the solution of the 
Reynolds equation for the pressure field. The second case includes the effects 
of fluid inertia at the film lands and a "ram" pressure effect at the pad 
leading egde (Kp=0.25). Referring to Figure 6 of Taniguchi et al. paper (1988), 
the present inertial flow model correlates better with the test results for 
pressure and film thickness. It is noted here that the numerical predictions 
without fluid inertia are virtually identical to those of Taniguchi et al. 
Recently, Fillon et al. (1991) performed similar comparisons with Taniguchi's 
test results and showed better agreement when accounting for the effects of 
elastic distortion due to thermal growth and pressure loading. 

The film temperatures presented in Figure 9 are low compared with the 
temperature measurements at the pad bearing surface. The error is as large as 
10°C for the loaded pads. Note that the present model predicts a film bulk-flow 
temperature , i.e. an average temperature across the film thickness, which is 
certainly smaller than the measured temperatures on the bearing pads or 
predicted across the film thickness in a 3-dimensional model. Nonetheless, the 
results of the present analysis are shown to predict very well the static 
equilibrium journal position and provide an accurate representation of the film 
thickness variation with a representative distribution of the film bulk-flow 
temperature. 

Someya's journal bearing databook (1988) is a standard reference for the 
design of fluid film bearings. The book details the theoretical (numerical 
calculations) static and dynamic force performance characteristics for a variety 
of fluid film bearing geometries, and includes test results for a number of 
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relevant bearing configurations. Test bearing No. 11 of Someya's scholar 
reference, namely a tilting pad bearing with five pads, has been chosen to 
validate the numerical predictions from the present analysis. Table 2 shows the 
geometric and operating characteristics of the test bearing article. Someya 
presents the test results in a dimensionless form as a function of a Sommerfled 
number defined as S^SJ-L-D (R/Cy [ 2itWJ where the effective viscosity is 
assumed to have the value corresponding to the measured fluid outlet 
temperature. Calculations using the present model have been performed for two 
flow conditions, a) an isothermal flow model without thermal energy transport 
and using the effective viscosity provided by Someya; and, b) an adiabatic flow 
model with null heat flow through the bearing and shaft surfaces, a heat 
carry-over mixing coefficient equal to 0.50, and an oil viscosity vs. 
temperature relation based on an exponential formulae. The test bearing operates 
on the laminar flow regime with no appreciable effects of fluid inertia at film 
lands or pad leading edges. The isothermal model was carried out to compare 
results with the analytical predictions presented by Someya (1988). 

Figure 10 shows the measured and predicted equilibrium eccentricity 
versus the bearing Sommerfeld number. The largest Sommerfeld number (S=0.106) 
denotes test conditions of low load (1.2 kN) and high shaft speed (6,120 rpm), 
while the lowest Sommerfeld number represents the largest load (9.48kN) and 
lowest journal speed (2,982 rpm). The isothermal results are virtually identical 
to those reported by Someya (1988) in his handbook. Note that the adiabatic flow 
model provides better agreement with the test data for the low speed, high load 
condition (small S) where thermal effects are expected to be of importance. 

The experimental synchronous (direct) stiffness and damping force 
coefficients obtained by Someya are given in Figures 11 through 14. The 
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analytical values obtained from the isothermal flow model are similar to those 
from the theoretical model of Someya and not reproduced for brevity. The figures 
show a different scaling from those depicted on the databook to emphasize the 
differences between tests and theory. Figures 11 and 12 show the dimensionless 
direct stiffness coefficients (K^i 
respectively. The adiabatic flow model predicts better the stiffness 
coefficients in particular at low speeds and high loads (low Sommerfled 
numbers). The large differences between the isothermal model and the adiabatic 
model at low S's are due to the heating of the fluid film at large loads which 
brings a temperature rise as large as 100°C (calculated) on the film lands. The 
experimental values for the direct stiffness are very low compared with the 
theoretical results. The large differences are explained in terms of both 
thermal and elastic deformations since the largest predicted pressure is 20.6 
MPa (2,987 psi). 

Figures 13 and 14 present the dimensionless damping coefficients (C^y, 
Cjy) versus the effective Sommerfeld number. The predicted damping coefficients 
agTee very well with the test data, while the other direct coefficient 
shows poor correlation at the lowest values of Sommerfeld number (large load, 
low speed). The theoretical results for the adiabatic flow model correlate 
better and at least do not increase as the Sommerfeld number decreases. 

It is well known that most analyses overpredict damping coefficients 
when compared to test results in tilting-pad bearings. The rationale for the 
large discrepancies is not yet clear although thermoelastic pad deformations and 
support structural effects are suspect (Nicholas and Barrett, 1986, Kim and 
Palazzolo, 1993). The comparisons shown in the figures demonstrate the current 
model to agree well with test results for low loads and high speeds conditions. 


K^y) versus the Sommerfeld number, 
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Thermoelastic pad deformations could be a major source of discrepancy for large 
loads. Note that the predictions from the present model are representative of 
other models available in the literature. 

Predictions obtained with the present computational program are 
identical to those presented by Armentrout and Paquette (1993) and Chen (1994) 
for flexure-pad bearings operating at high speeds and light viscous fluids under 
laminar flow conditions. Typical results show, at the design load condition and 
as as the pad rotational stiffness increases, the variation of force 
coefficients from a tilting-pad regime passing through a transition zone and 
leading to a fixed-geometry regime. Designers typically select the value of 
rotational stiffness which shows the least value of the difference in cross - 
coupled stiffness coefficients (K^ K^) just before the transition regime 
starts. Then a finite element program is used to design the geometry of the 
flexure pad web which renders the required pad rotational stiffness as well as 
pad radial stiffness when the case merits. In its simplest form, the flexure pad 
tational and radial stiffness are given by the following elementary relations: 

V E *V h w l N * m /radJ; 1^ = A w -E/h w [N/mJ (56) 

where E is the elastic modulus of the web material, h is the web height, A =w-L 

W W 

is the web cross-sectional area (of width v and axial length L) , and I =L*w*Vl2 

w 

is the area moment of inertia of the web. The expressions detailed do not 
account for the geometric stiffness due to the load applied on the pad. Note 

that the ratio of the radial stiffness to rotational stiffness (1^/K ) is equal 
to 12/w 2 . 

The following discussion pertains to the results obtained for a flexure 
pad hybrid bearing for application in a liquid oxygen (LOx) turbopump. The 
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geometry and operating conditions for the six pad bearing are given in Table 3. 
Each pad has a rectangular recess and an orifice fluid supply line machined 
across the web supporting the flexure pad. The recess to pad ratio of areas is 
equal to 0.236 as recommended for a LOx application (Butner and Murphy, 1986). 
The values of supply and discharge pressure and rotational speed correspond to 
that on an ALS bearing turbopump configuration. The large values of axial flow, 
circumferential, and squeeze film flow Reynolds numbers show a turbulent flow 
application with large fluid inertial effects. Calculations for both isothermal 
conditions and adiabatic heat flow with full pad leading edge thermal mixing 
were performed simultaneously. The resuls between both models differ little 
expect for an increment in fluid temperature at the pad sides of 12°K and not 
sufficient to affect considerably the bearing performance. Henceforth, in the 
results that follow only the isothermal results are discussed in detail. 

At the journal concentric position (e=0), Figure 15 shows the calculated 
whirl frequency ratio (WFR) and recess pressure ratio for increasing values of 
the pad rotational stiffness. The lowest value of represents an ideal tilting 
pad geometry with null restraining moments, while the largest rotational 
stiffness effectively represents a rigid (fixed) pad bearing configuration. The 
effect of pad angular stiffness on the stability indicator (WFR) is clear and 
shows the advantage of a moving pad geometry for this high speed application. As 
the pad rotational stiffness increases, the recess pressure ratio decreases 
while the flow rate increases as shown in Figure 16. The (X) marked on this 
figure represents the flow rate for a conventional fixed HJB with film lands of 
360 extent. Figure 17 shows the pad rotation angle and the pad minimum film 
thickness (occurring at the pad trailing edge) versus the pad rotational 
stiffness. As this stiffness increases it is clear that the pad is unable to 
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m ° Ve and he " Ce the n,ini " un ' £il " thickness is that of the fixed pad. The largest 
rotations are expected in the tilting pad regime with a reduction in film 

thickness of 40*. The drag torque shown in Figure 18 is shown to decrease with 
the pad rotational stiffness since the effective film thickness is essentially 
larger. This result then shows that the flexure pad bearing has lower drag power 
than the conventional tilting pad bearing. 

Figure 19 shows the synchronous force coefficients, i.e. evaluated at 
<*=B. At the concentric position, the direct coefficients are equal ) 

while the cross-coupled coefficients differ in sign (K^-K^). Note that the 
direct stiffness is about an order of magnitude larger than the cross-coupled 
stiffness coefficients for the fixed pad condition. The calculations show the 
paramount effect of the pad rotational stiffness coefficient on the rise of the 
cross-coupled coefficients. Figures 20 to 23 show the effect of frequency 
excitation on the stiffness (K^, K^), and damping (C^, C^) coefficients, 
respectively. The results denote calculations at frequency ratios /.«/ fi equal to 
0, 0.5, 1 and 2 times the rotational speed. The direct stiffness decreases 
with the frequency ratio / mainly due to fluid inertia effects, with a maximum 
direct stiffness at a pad rotational stiffness equal to 10,000 N-m/rad. The 
excitation frequency appears not to affect the cross-coupled coefficients 
showing then low values of cross-coupled inertia coefficients. 

Next, a pad rotational stiffness (KJof 10 KNm/rad was selected to 
perform calculations for the load capacity of the bearing at the rated operating 
condition. The load (V) applied to the bearing increased to a maximum of 12 kN 
(2,700 lbs) denoting a specific load (W/LD) equal to 3.5 MPa (507 psi). Two 
cases were considered, (A) load applied towards the center of a bearing pad - 
LOP, (B) load applied at 30° from the vertical line, i.e. between pads - LBP. 
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Figure 24 shows the calculated journal eccentricity and the minimum file 
thickness occurring on the bottom pad to have a very linear relationship vith 
pplied load. Note that the smallest film thickness is just 402 of the 
nominal bearing clearance for the largest load applied. The selected flexure pad 
rotational stiffness has a whirl frequency ratio equal to 0.25. This pad 
Stiffness value although not low enough to elliminate hydrodynamic inestability 
constitutes a major advancement over the rigid bearing pad configuration. 

The flexure pad geometry to provide the selected value of pad rotational 
stiffness can be obtained using the formulae given by equation (56). Using a 
conventional steel with an elastic modulus of 20.7E10 N/„ 2 (30E6 psi), a web 

ght ( V ° £ 1 cm and thickness (w) of 5.4 mm renders a value K^IOKN- m/rad. 

The web thickness is enough to allow for the orifice supply line of 2.581 mm 
diameter. The pad radial stiffness is equal to 4,135 MN/m which is an order of 
magnitude larger than the fluid film direct stiffness K^, and hence it will 
not affect greatly the bearing performance. (Maximum elastic web deformations of 
3um are expected for the largest load applied of 12KN). 

Figures 25 and 26 show the synchronous stiffness and damping coefficients 
versus the applied load, respectively. Note that the stiffness coefficients are 
virtually constant for loads to 8kN (1800 lbs). The cross-coupled stiffnesses 
for the case of load between pads (B) show the largest deviations at the largest 
load applied, while the direct damping coefficients decrease slightly with the 
applied load. The whirl frequency ratio decreases monotonically from 0.25 to 
0.22 as the load increases. Figure 27 shows the calculated (dimensionless) 
centerline pressure and film thickness for a load of 12 kN. The largest film 
pressures along with minimum film occur at pad 4 where the journal load is 
applied. The results show a small "ram" pressure at the pads inlet and 
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hydrodynamic pressure increases in the downstream portion of the bearing 
recesses. Note also the large inertial pressure drops at the recess edges. 

The analysis also includes a pad rotational damping and set to zero in 
all calculations presented. Flexure-pad bearings due to its automated wire-EDM 
construction offer the advantage of accurate specification of the gap (g) 
between the pad inersurface and the bearing housing. A simple one-dimensional 
analysis shows that the pad rotational damping is approximately equal to: 

C r = (l/60)-p-L-(9 p -nD/180) 5 (l/g) 3 [N-m-sec/rad] (57) 

The rotational damping for a gap (g) equal to 101.6pm (0.004 mils) is equal to 
0.05187 N-m*sec/rad and not large enough to affect the dynamics of the LOx 
flexure pad bearing. However, some commercial applications handling high 
viscosity lubricants may offer large values of rotational damping. It this 
damping is large enough it will cause the flexure pad to become overdamped and 
in fact, it may lock the pad and degrade the dynamic performance of the whole 
bearing. This condition is very interesting since the flexure-pad bearing will 
still ‘behave well under static load conditions, i.e. it will show a journal 
displacement in the same direction as the applied load, but due to excess in 
rotational damping, it can bring hydrodynamic stability at sufficiently large 
speeds and whirl frequencies. 

The calculations performed show that a flexure-pad bearing keeps 
virtually all the benefits of a hydrostatic pad application (large direct 
stiffness and damping coefficients) while offering an accurate control on the 
cross-coupled stiffness coefficients and reduction in the whirl frequency ratio. 
The results provide ample evidence on the effect of pad rotational stiffness on 
the performance of a hybrid - flexure pad bearing. 
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CONCLUSIONS 


This reports details the theoretical analysis for calculation of the 
dynamic force performance of tilting-pad and flexure-pad hybrid (combination 
hydrostatic and hydrodynamic) bearings operating in the turbulent flow regime. 
These types of bearings appear to offer an alternative to improve (and 
eliminate) the limited stability characteristics of conventional fixed pad 
geometry bearings. In particular, a flexure pad bearing is best suited to meet 
this objective since it offers many advantages when compared to a conventional 
tilting pad bearing. These advantages encompass acuracy of manufacturing without 
added difficulties for assembly and calibration, adequate control of bearing 
preload and pad structural stiffnesses (radial and rotational), and absence of 
wear between shoe and pivot support. Flexure pad bearings lend themselves to a 
hydrostatic application since the ports for the pressurized fluid can be easily 
machined on the pad supporting web without added mechanical complexities and at 
a reduced cost. 

Bulk-flow equations of mass conservation, momentum and energy transport 
describe the motion of a variable properties fluid within the thin film lands 
of a bearing. A mass conservation equation at the recess volumes with a simple 
formulation for the recess pressure and temperature fields are also considered. 
Fluid inertia effects, temporal and advective, are fully accounted for flows 
with large pressure differentials or large surface speeds. Zeroth- and 
first-order flow equations are numerically solved to determine the flow field 
and calculate the pad fluid film forces and force coefficients. An elegant 
method to determine the journal equilibrium position and the reduced force 
coefficients for tilting-pad bearings is introduced. 
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Numerical results from the present development have been compared to 
similar analyses for laminar flow tilting-pad bearings. Predictions from the 
model also correlate well with experimental bearing loads and dynamic force 
coefficients. The comparisons insure the validity of the model implemented and 
also widen the range of practical applications for the computational program 

developed. 

Predictions for a flexure-pad hybrid bearing handling liquid oxygen at 
operating conditions typical of a turbopump applications demonstrate the ability 
of this type of bearing to improve the hydrodynamic stability whirl ratio, 
without degradation on the static load performance or reduction in direct 
stiffness and damping coefficients. The test example shows that hybrid bearings 
satisfy in excess the load requirements typical of present and future cryogenic 

turbopumps . 

The analysis still needs to be improved to account for pad flexibility 
in the radial direction. Industrial applications with large loads applications 
may also demand a thermoelastic model for accurate determination of pad surface 
temperature and minimum film thickness. These considerations may not be of great 
importance in a cryogenic turbopump application. 
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Table 1. Geoaetry and Operating Condition of Taniguchi et al. Bearing (1990) 

Number of pads = 4 

Diameter (D) = 479 mm 

Length (L) = 300 mm 

Clearance (C p ) = 612 pm, zero preload 

Pad Arc: 80°, offset=0.50 (centrally pivoted) 
pivots at 45°, 135°, 225°, 315° 

Pad thickness=121 mm 

Lubricant: Oil ISO VG 32 with p= 860 kg/m 3 
viscosity y=0.0277 Pa-s at 40°C 
oil viscosity-temperature coefficient 0.034 l/°c 

Supply conditions: Ps=0.1 MPa at Ts=40°C 

flow rate 520 lit/min 

Operating speed: lkrpm to 4krpm 

Load between pads 2 & 3: 60kN to 300 kN. 

Circumferential flow Reynolds number (pQRC/y) = 1,430 
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Table 2. Geometry of Test Bearing No. 11 of Someya (1988), pp. 227-229. 

Number of pads = 5 

Diameter (D) = 76.4 mm 

Length (L) = 31.6 mm 

Clearance (Cp) = 86 ym, zero preload 

Pad Arc: 58°, offset=0.50 (centrally pivoted) 
pivots at 36°, 108°, 180°, 252°. 324° 

Lubricant: #140 turbine Oil with p=874 kg/m 3 

viscosity v= 61.8 mm 2 /sec at 37.8°C 
35.4 " 50.0°C 

7.78 ' ' 98.9°C 

Supply conditions: Ps=0.1 MPa 


Inlet Temperature: varies between 38.5°C to 41.3°C 
Operating speed: varies between 2,980 rpm to 6,240 rpm 

Load V on bottom pad (#3): 1.2 kN to 9.48 kN. 


Sommerfeld Number 



S L D 


V 



2 


Dimensionless 


Refer to Someya' s 


stiffness coefficients k„ = K„ (C/V) 
damping "» c _ = c (C-Q/W) 

Databook for a detailed description of test 


conditions 


- o - 
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Table 3. Flexure-pad Liquid Oxygen Hybrid Bearing. 

Geometry and operatinj conditions 

Diameter D=92.71 mm (3.65 in) 

Length L=37.08 mm (1.46 in) 

Clearance C = 76.2 ym (0.003 in), zero preload 
Recess depth H R = 228ym (0.009 in) 

6 pads of length 50° and 1 recess per pad 
pad offset=0.50 

pad pivots at 0° , 60°, 120° , 180° , 240°, 300° 

leading edge ram pressure coefficient K p =0.32 



recess arc length 23.5° and axial length 1= 19mm (0.75 in) 


orifice diameter d Q =2.58 mm, C d =1.0 

edge loss coefficients =0.0 

swirl ratio ot=0.50 W ^ 


Pad rotational stiffness (K f ) varies 

damping (C r )=0 N-m-s/rad 

Inertia : 1.017E-4 kg-m 2 , mass : 0.283 kg for thickness 0.02 m 
Operating speed: 25krpm 

Pressure supply P = 26.71 MPa (3874 psi) 

discharge P® = 8.81 MPa (1278 psi), AP= 2600 psi 

Liquid oxygen at inlet temperature T =90°K (162°R) 

s 

at < P s-V : densit y P s = 1 > 192 kg/m 3 , viscosity u s =0. 2459E-03 Pa-s, 

Reynolds numbers at concentric position. 

Circumferential flow Re =(p/y) G-R-C= 44,832 

Axial flow Re^=(p/y) a V-C = 30,487 - 30,801 

Squeeze film Re =(p/y) w-C = 73.52 (<*=Q) 

s s 

- o - 
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Figure 1. Geometry of flexure pad hybrid bearing 
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preload rp_ = r P /Cp ; 
rp=Cp-Cm 

r&_=0: cylindrical pad, 
m_=l ; journal - pad contact 


t > 


t R+Cp 


pivot 


R+Cm 








. . ... ........ ■ • /- I 



Cp: pad clearance j 
Cm: assembled clearance 


film thickness 


H=Cp+ex cos(0)+ev Sin(6)-r p cos(0-0 P )-R 5 sin(0— 0 P ) 


Figure 2. Geometry and nomenclature for 

tilting-pad 
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Figure 3- Conceptual description of pressure rise 

and drop at recess edge of a hydrostatic 
pad bearing, and pressure ram effect at 
leading edge of bearing pad.. 
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ISOTH=5, Isothermal journal & Radial heat flow through bearing 
ISOTH=4, Adiabatic journal & Radial heat flow through bearing 


Qradial 



her = fiB (Kb/R)/{ Kb/R + hs In(RouVR)} 
where Jib is the fluid- film to bearing surface 
convection heat transfer coefficient 


Figure 4b. Radial heat flow on bearing shell 
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Janiguchi 4pad bearing, 3Krpm 
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Load, Wx (KN) 

Figure 5. Equilibrium eccentricity vs. applied load 

Taniguchi et al. tilt pad bearing. Speed 3krpm 




Taniguchi 4pad bearing, 3Krpm 
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Load, Wx (KN) 

Figure 6. Minimum film thickness (|im) vs. applied load 

Taniguchi et al. tilt pad bearing. Speed 3krpm 
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Figure 7. Centerline film pressure (MPa) on loaded pads 

Taniguchi et al. tilt pad bearing. Speed 3krpm, Load=l80kN 



Taniguchi 4pad bearing, 3Krpm, 180kN - LBP 



£ $ Z spBd uo (wn) sseu>jo!L|i tuiy 


67 


Circumferential cooordinate (degs) 

Figure 8. Film thickness (jim) on loaded pads. 

Taniguchi et al. tilt pad bearing. Speed 3krpm, Load=1 80kN 



Janiguchi 4pad bearing, 3krpm, 180kN - LBP, Ts=40C 
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Tilt pad bearing, Someya HB, LOP 
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Figure 10. Equilibrium eccentricity vs. Sommerfeld Number 

TEST Measurements, isothermal and ADIAbatic solutions. 
Someya s databook 4 shoe tilt pad bearing - Load on pad 



Tilt pad bearing, Someya HB, LOP 
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Figure 11. Direct Stiffness Kxx[C/W] vs. Sommerfeld Number 

TEST Measurements, isothermal and ADIAbatic solutions. 
Someya's databook 4 shoe tilt pad bearing - Load on pad 



Tilt pad bearing, Someya HB, LOP 
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Figure 12. Direct Stiffness Kyy[C/W] vs. Sommerfeld Number 

TEST Measurements, isothermal and ADIAbatic solutions. 
Someya’s databook 4 shoe tilt pad bearing - Load on pad 



Tilt pad bearing, Someya HB, LOP 
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Figure 13. Direct Damping Cxx[C£2/W] vs. Sommerfeld Number 

TEST Measurements, isothermal and ADIAbatic solutions 
Someya's databook 4 shoe tilt pad bearing - Load on pad 




Tilt pad bearing, Someya HB, LOP 
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Figure 14. Direct Damping Cyy[CO/W] vs. Sommerfeld Number 

TEST Measurements, isothermal and ADIAbatic solutions. 
Someya s databook 4 shoe tilt pad bearing - Load on pad 




LOX 6 flexure pad HJB, 25Kcpm, Ps-Pa=17.9MPa, T=90K 
concentric e=0, do=2.581mm 
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Figure 15. Whirl frequency ratio and recess pressure ratio vs. 

pad rotational stiffness. LOx flexure-pad hybrid bearing 
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LO c X on 6 cl e nW=!)^o J =llife m ' p s-Pa=17.9MPa, T=90K 
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Figure 17. Pad rotation and minimum film thickness vs. 

rotational stiffness. LOx flGxure— pad hybrid bGarinQ 
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Pad Rotational Stiffness [N.m/rad] 

Figure 18 . Drag Torque vs. pad rotational stiffness. 

LOx flexure-pad hybrid bearing - concentric operation 
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Pad Rotational Stiffness [N. m/rad] 

Figure 20. Effect of frequency on direct stiffness Kxx=Kyy 
LO x flexure-pad hybrid bearing - concentric operation 



concentric e=0, do=2.581mm 
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Pad Rotational Stiffness [N. m/rad] 

Figure 21. Effect of frequency on cross stiffness Kxy=-Kyx 
LO x flexure-pad hybrid bearing - concentric operation 



LOX 6 flexure pad HJB, 25Kcpm, Ps-Pa=17.9MPa, T=90K 
concentric e=0, do=2.581 mm 
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Pad Rotational Stiffness [N. m/rad] 

Figure 22. Effect of frequency on direct damping Cxx=Cyy 
LOx flexure-pad hybrid bearing - concentric operation 




LOX 6 flexure pad HJB, 25Kcpm, Ps-Pa=17.9MPa, T=90K 
concentric e=0, do=2.581mm 
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Pad Rotational Stiffness [N. m/rad] 

Figure 23. Effect of frequency on cross damping Cxy=-Cyx 

LOx flexure-pad hybrid bearing - concentric operation 



LOX 6 flexure pad HJB, 25Kcpm, Ps-Pa=17.9MPa, T=90K 
A. Load on pad 4, B: Load between pads 4&5 
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LOX 6 flexure pad HJB, 25Kcpm, Ps-Pa=17.9MPa, T=90K 
A. Load on pad 4, B: Load between pads 4 & 5 
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Figure 27. 
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ABSTRACT 

An analysis for the prediction of the transient response of a rigid rotor 
supported in turbulent flow, externally pressurized fluid film bearings is 
presented. The rotor-bearing system nonlinear equations of motion are 
numerically solved with local linearization at each integration time. Fluid film 
bearing reaction forces are calculated from the numerical solution of the 
unsteady bulk-flow equations on the film lands. Analytical perturbations 
determine linear bulk-flow equations and their solutions allow estimation of 
rotordynamic force coefficients. The current numerical model is largely based on 
prior efficient and accurate computational programs for the prediction of the 
static and rotordynamic force response characteristics of turbulent flow fluid 
film bearings. 

Examples for the transient response of rigid rotors supported on plain 
journal bearings, annular pressure seals and hydrostatic bearings under a 
variety of external loads are discussed in detail. Transient responses based on 
an approximate model which uses rotordynamic force coefficients at an 
equilibrium position are also presented. The full nonlinear model demands small 
time steps and is computationally very intensive. For externally pressurized 
bearings, the linear transient model provides good results in a fraction of the 
time required for the nonlinear model to compute the rotor-bearing dynamic 
forced response. 


For a lucid discussion on the importance and advantages of fluid film bearings 
in cryogenic turbomachinery see the Research Progress Report, Year I, to the 
present project (San Andres, 1993) 


Refer to hydrotran User's Manual and Tutorial for information and instructions 
to run the program. 
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NOMENCLATURE 

A 


2 

Cd n d Q / 4 . Equivalent orifice area [m 2 ] 


R ecess area [m 2 ]. 

Rotor accelerations [m/s 2 ]. 

Radial clearance function [ m ] . 
p Fluid specific heat [J/kg-°K]. 

C XX’ C XY’ C YX’ C YY Dam Pi n g force (linear) coefficients [Ns/m] 


A X ,A Y 

C 

C 


D 

d 


J,B 


F F 
X’ Y 


Orifice discharge coefficient 
2-R . Bearing diameter [m]. 

Orifice diameter [m] 

e u 

; a =0.001375 
b =500,000; c M =10,000 
eJ]=l/3.00 " 

Turbulent flow friction factors at journal and bearing surfaces. 
Bearing fluid film forces along {X,Y} axes [N]. 


1 + 

C 1±’Z + b M 1 

M 





1 H r j,b J 

■ 


H 

H 


h B’ h J 


K XX’ K XY' K YX’ K YY 


K 

eq 

L,1 

M 


m xx ,m xy ,m yx ,m yy 

M 

N 

recess 
P, Pr 


C + X(t) cos ( 0) + Y(t) sin(0) . Film thickness [m] . 

Recess depth [ m ] . 

heat transfer coefficients to bearing and journal [V/m^°K] 
Bearing force stiffness coefficients [N/m] 

Equivalent linear bearing stiffness [N/m] 

Bearing axial length, recess axial length [ m ] . 

Point rotor mass [kg]. 

Bearing inertia force coefficients [kg]. 

VV p r * 0rifice circumferential velocity Mach number. 
Number of recesses on bearing. 

Fluid pressure, recess pressure [N/m 2 ]. 

Pressures just before and after recess edge [N/m 2 ]. 
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P 0' P L 


Re 

Re 



r 


J’ 


r 


B 


S. 

1 




T,T 


U ,U 
x T y 


V 

s 

V V Y. 

w x (t),w Y (t) 

X(t),Y(t) 


Z 


&r 

n 

e 


0 

r 

0 


Discharge pressures on left and right sides of bearing[N/m 2 ] . 
External supply pressure [N/m 2 ]. 

heat flow to bearing and journal surfaces [U/m 2 ]. 

Bearing radius [m]. 

(p-G-C-R/p)*. Nominal circumferential flow Reynolds number. 

2 

(p wC /vO*. Nominal Squeeze film Reynolds number 
(p/m)hJ[(U 9 -G-R) 2 + U z 2 ]; (p/y)Hl[U e 2 + U z 2 J 

Flow Reynolds numbers relative to journal and bearing surfaces. 
Roughness depths at journal and bearing surfaces [m]. 

0'Atj. Extended time step [ s ] • 

time coordinate [s], initial time [sj. 

Temperature, supply temperature [°KJ. 

Bulk-flow velocities in circ.(x) and axial (y) directions [m/s]. 
(H r +H)A r +V s Total recess volume [m^]. 

Volume of orifice supply line [m 3 ]. 

Journal (rotor) velocity components [m/s]. 

External loads to rotor-bearing [N], 

journal (rotor) displacements in inertial coordinates [m]. 

T 

[X,Y] . Vector of journal displacements [m]. 

(l/p)( 3p/9P) . Liquid compressibility coefficient [m 2 /N]. 

-(l/p)( 3p/3T) . Liquid volumetric expansion coefficient [ 1 / °K ] . 
^(H r +H). Ratio of land film thickness to recess depth. 
x/R. Circumferential or angular coordinate, 
rth-recess angular length [rad], 
time integration parameter, typ=1.42. 
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K =K 
Z X 


V k b 


P»U 


k xd 


S2, « 

<\> 

T 

T 

n 


r 

At 


AZ. 

1 

AV Zi 

M Zi 


2 (K J +K B ) - Tur bulence shear factors in (y,x) flow directions. 

f j' R J’ f B " R B ’ turbulent shear parameters at journal and 
bearing surfaces. 


Re 


0.681 


7.753 


Turbulent shear flow parameter at recess. 


Fluid density [Kg/m 3 ], viscosity [Ns/m 2 ]. 

® r ^«Jsr^rSrM^r* £fici ““ in circu " £ — 

direction. recess_ed ® e entrance loss coefficients in axial 

Rotational speed of journal, excitation or whirl frequency [1/s] 
whirl frequency ratio. 


t t^. Local time coordinate [ s ] - 

2n/ V Natural P er iod of rotor-bearing system [s]. 

Q-R. Journal surface velocity [m/s], 
recess boundary with outward normal n. 

Time step [s] 

Z i+l~ Z i * Vector of incremental displacements [m]. 

V Zi+l~ V Zi- Vector of incremental journal velocities (m/s] 
A Zi + l“ A zr Vector of incremental journal accelerations [m/s 2 ] 


Subscripts refer to: 




r,e 




B, J 


° f l0Cal circumfer ential and axial 
plane of bearing. 

bearing recesses and edges (entrance). 


coordinates 


in 


refer to discrete times t. and t 

J i* 

refer to bearing and journal surfaces. 


90 



LIST OF TABLES 


Table 1 


Table 2 


Table 3. 


Table 4. 


Description of rotor supported on plain cylindrical 
earings. Transient response due to applied static 
weight. Data from Childs et al. (1977). 

Description of rotor supported on plain cylindrical 
earings . Transient response due to sudden unbalance 
load* Data from Childs et al. (1977). 


Description of rotor supported on plain cylindrical 
bearings. Transient response due to impulse load. 
Data from Tichy and Bou-Said (1991) . 


Description of rotor supported on damper seal 
bearings. Seal data from San Andres and Yang (1994) . 

Description of a rotor supported on hydrostatic 
bearings. Bearing data from Childs and Hale (1994) . 


Table 5. 



LIST OF FIGURES 


Figure 1. Rigid rotor supported on fluid film bearings. 

Figure 2. Geometry of a Hydrostatic Bearing. 

Figure 3. Conceptual description of pressure rise drop 
et recess edge of a hydrostatic pad bearing, and 
pressure ram effect at leading edge of bearing pad. 

Figure 4. Linear acceleration in numerical integration method. 

Figure 5. Transient response of Rigid Rotor on Cylindrical 

Journal Bearings. Sudden step load Wx«Mg=8,000 [N] . 

Figure 6. Fluid film forces and inertia forces for transient 
response. Sudden step load Wx=Mg=8, 000 [N] . Wy=0. 

Figure 7 . Transient response of Rigid Rotor on Cylindrical 
Journal Bearings. Sudden unbalance load 
Wx=Mg+M u to 2 cos (cot), Wy=M u CD 2 sin (cot); u=25.4um 

Figure 8. Journal displacements of Rigid Rotor on Cylindrical 
Journal Bearings. Sudden unbalance load 
Wx=Mg+M u CO 2 cos (cot), Wy=M u CO 2 sin(COt); u=25.4um 

Figure 9. S— S fluid film forces for transient response. 

Sudden unbalance load Wx=Mg+M u co 2 cos (cot), 

Wy=M u CO 2 sin(cot); u=25 . 4um 

Figure 10. Effect of fluid inertia on transient response. 

Impulsive load Wx=M. g=2kN, Wy=-20kN over 2.77 msec. 

Figure 11. Journal displacements with and without fluid inertia. 

Impulsive load Wx=M.g=2kN, Wy=-20kN over 2.77 msec. 

Figure 12. Bearing forces with and without fluid inertia. 

Impulsive load Wx=M. g=2kN, Wy=-20kN over 2.77 msec. 

Figure 13. Rotor inertia forces with and without fluid inertia. 

Impulsive load Wx=M. g=2kN, Wy=-20kN over 2.77 msec. 

Figure 14. Transient seal locus due to ramp and step loads. 

Static load Wx=10,000 N, Wy*0 N. 

Figure 15 . Transient LOX seal eccentricity due to step r amp 
loads. Static load Wx=10,000 N, Wy=0 N. 

92 



LIST OF FIGURES (continued) 


Figure 16. Transient bearing force Fx due to step and ramp 
loads* Static load Wx=10,000 N, WyaO N. 

Figure 17 . Transient seal center locus due to inclusive load. 
Wx=M.g=490 N, Wysio kN (impulse), to=0.57 ms. 

Figure 18. Transient seal force Fy due to impulse load Wy 
Wx=M. g=490 N, Wy=10 kN (impulse), to=0.57 ms. 

Figure 19. Transient seal eccentricity due to periodic load. 
Load 5 kN synchronous (433.3 Hz). 

Figure 20. Transient bearing force due to periodic load. 

Load 5 kN synchronous (433.3 Hz). 

Figure 21. Transient seal center locus due to periodic load. 
Load 5 kN synchronous (433.3 Hz). 

Figure 22. Transient HOB journal center locus due to ramp loads 
Load Wx=5 (t/T*) kN, T*=2 msec, 0.2 msec, Wy=0 N. 

Figure 23.HJB journal eccentricity for ramp loads. 

Load Wx=5 (t/T*) kN, T*=2 msec, 0.2 msec, Wy=0 N. 

Figure 24. Transient bearing force Fx due to ramp loads. 

Load Wx=5 (t/T*) kN, T*=2 msec, 0.2 msec, Wy=0 N. 

Figure. 25. Transient journal locus due to impulsive load. 

Wx=M.g=490 N, Wy=5 exp(-t/to) kN, to=0.5 ms. 

Figure 26. Transient journal center displacements for impulse 
load. Wx=M. g=490 N, Wy=5 exp(-t/to) kN, to=0.5 ms. 

Figure 27. Transient HJB force Fy due to impulsive load Wy. 

Wx*M. g-490 N, Wy=5 exp(-t/to) kN, to=0.5 ms. 

Figure 2 8. HJB journal eccentricity for unbalance load. 

Wx=Mg+M . U . to 2 cos ( cot ) , Wy=M . U . CO 2 sin(tot), U varies. 

Figure 29. HJB steady state orbit for umbalance load. 

Wx=Mg+M . U . co 2 cos ( cot ) , Wy«M.U.C 0 2 sin(cot), U varies. 

Figure 30. Initial X-forces for unbalance U=40um, Fumb=19.71 KN. 

Figure 31. Steady-state X-forces for unbalance Us40um, 

Fumb=19 .71 kN. 



INTRODUCTION 


The transient forced response of rotor systems supported on fluid film 
bearings is of importance to study (or to discover) the effect of bearing 
nonlinearities brought by extreme loading conditions or stringent operating 
conditions. There is a general belief that classical rotordynamic studies based 
on linearized bearing force coefficients are unable to predict system time 
responses for realistic levels of external loads. Some types of fluid film 
bearing elements such as hydrodynamic bearings, foil bearings and squeeze film 
dampers are known to be highly nonlinear, and therefore, nonlinear transient 
analyses are necessary to establish the ranges for safe operation of these rotor 
bearing systems. Transient studies including accurate finite length fluid film 
bearing models can be performed in reasonable amounts of time thanks to the 
current high speed computing facilities. Furthermore, novel and sound analytical 
techniques for handling highly nonlinear systems are at the researcher's 
disposal to study aperiodic responses, jump-phenomena, bifurcation and even 
"chaotic" responses in certain regimes of operation (Ehrich, 1992). 

The development of an all fluid film bearing technology for cryogenic 
turbopumps has advanced steadily over the past years. Low cost, reliable fluid 
film bearing elements are needed to provide maximum operating life with optimum 
and predictable rotordynamic characteristics. San Andres (1990 - 1993) provides 
the analysis and computational programs for the calculation of static and 
rotordynamic force characteristics in turbulent flow, cryogenic liquid 
hydrostatic bearings and pressure seals. Extensive experimental tests performed 
at the high speed, high pressure, water - bearing test facility at Texas A&M 
University have validated the numerical predictions (Kurtin et al., 1993, 

Lindsey, 1993, Franchek et al., 1994). Further analytical and experimental 
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studies currently address the important issues of fluid compressibility with 
cryogenic liquids and the limited operating range due to hydrodynamic 
inestability in conventional bearing designs- The present investigation develops 
an analysis for computation of the transient dynamics and forced response of a 
rigid rotor supported on turbulent flow, fluid film bearings. The ultimate aim 
of the investigation is to provide a sound computational program for the 
prediction of synchronous unbalance response, rotor lift-off, and response to 
abnormal shock and maneuvering loads in cryogenic turbomachinery. 

LITERATURE REVIEW 

The literature related to the transient response of rotating machinery 
supported in fluid film bearings is quite extensive. From the numerous relevant 
references I have selected a few of those works which have directly shaped my 
understanding of the nonlinear rotor - fluid film bearing problem. Booker (1965) 
and later Childs et al. (1977) developed the mobility and impedances methods for 
solution of rotor-bearing problems involving transient responses. The impedance 
formulation applicable to rotordynamic studies expresses the bearing 
instantaneous reaction loads in terms of a single (squeeze) velocity vector in a 
particular rotating frame. This simple though fundamental idea allows fast 
evaluation of the fluid film bearing forces. The mobility and impedance methods 
are applicable to cylindrical bearings and were quite popular before the advent 
of the high speed, personal computer. The original methods were graphical in 
nature and provided the trained engineer with an opportunity to obtain quickly 
the dynamic response of rotors supported on fluid film bearings with a few 
calculations. 

Kirk and Gunter (1970-1976) and Barrett et al. (1977) pioneered the 
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numerical study of transient dynamics of simple rigid and flexible rotors 
supported on short length, cylindrical hydrodynamic bearings. The study 
confirmed the findings of the classical linear stability analysis of Lund (1965) 
and detailed the severity of a sudden blade loss in rotor - bearing systems. 
These early vorks have left a profound imprint on our current understanding of 
rotordynamics and shaped many later developments, most notably those relevant to 
the coupling of floating oil seals to rotor - bearing dynamics (Kirk, 1988, 

1994). These early fundamental references are rarely acknovledged today although 
their examples and methods have been replicated in numerous investigations 
published on the last decade. 

The analysis of the dynamics of rotors supported in squeeze film dampers 
( SFDs ) is perhaps the single topic which has received most attention in terms of 
its transient response. Current analyses of rotor-disk assemblies supported on 
SFDs are based on overly simplified analytical expressions for fluid film forces 
as derived from the short journal bearing model with the so called a-film 
cavitation assumption (Li and Taylor, 1987, Chen et al. 1993, Sundararajan et 
al. 1994). Transient responses based on direct numerical integrations, and 
lately computational solutions based on current nonlinear dynamics models with 
elegant mathematical methods, show the forced response of rotor-SFDs systems to 
highly nonlinear with extreme sensitivity to unbalance levels. Little effort 
has been placed on finding a combination of operating and design parameters 
Which will avoid the system undesirable response. On the other hand, the 
richness of the nonlinear (theoretical) behavior has been exploited to prove 
beyond practical limits the amazing accuracy of the analytical models. However, 
there is little physical evidence and practical experience attesting to the 
veracity of the theoretical and computational predictions. In the last four 
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years, 1990-94, over 20 journal and conference publications have presented 
extensive theoretical nonlinear transient studies of the problem. Only Zhao e, 
>1. (1994) provide fundamental experimental results .hich demonstrate that 
current theoretical-computational models fail to predict the observed 

performance of actual rotor-damper systems. It appears then that the u-fil„ 

short length SFD model used in these nnnUnm 

onlinear analyses has little relationship 

with reality. 

Choy et al. <1990a-b, 1991) present detailed numerical studies on 

nonnnear effects in plain journal bearings due to journal misalignment, levels 

of loading and thermal conditions. Host importantly, the authors develop 

ctions for the response of a simple rotor - journal bearing system 

due to impulse loadings and compare them with transfer functions obtained from 

lin “ bear1 " 8 £ °™ “•« 1 “*"ts. The numerical results shoe the extent 
of nonlrnear behavior to be of importance for journal operation a, lev 

eccentricities ehere nonlinearities shoe a mashed increase as the journal is 

splaced to a larger operating point. Adams et al.(1983, 1994) present 
transient studies of point mass rotors supported on cylindrical journal bearings 
and tiiting-pad bearings. The analyses address the phenomena of pad flutter a. 
subsynchronous frequencies and the appearance of chaotic phenomena under 
particular operating and loading condi, ions. Tichy and Bou-Said (1,91, report on 
the importance of fluid inertia effects on the initial transient response of 
simple rotors supported in hydrodynamic plain journal bearings and acted upon by 
impulsive loads. The method extends the application of the f Urn-averaged 
momentum equations (Elrod et al., 1,83, to analyse unsteady floes in fluid film 
gs. This type of average-inertia model has been used extensively to 
address fluid inertia effects in squeese film dampers (Sseri et al. , 1,83). 
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Paranjpe and co-workers have advanced steadily the thermohydrodynamic 
analysis of dynamically loaded journal bearings for automotive applications. The 
latest publication (Paranjpe et al. , 1994) describes the complexity of their 
current model and provides a lucid explanation on the time scales for thermal 
effects in a fluid film bearing and its bounding solids. Desbordes et al. 

(1994), and Gadangi et al. (1994a, b) include the effects of film temperature and 
pivot flexibility on the transient response to suddenly applied unbalance loads 
in rotors supported on tilting-pad bearings. The authors bring to attention a 
rotor-bearing highly nonlinear behavior for large levels of unbalance, and where 
a linear rotordynamic model fails to predict the system transient response. 

For turbulent flow applications in annular pressure seals, Nordmann and 
Dietzen (1988), Tam et al. (1988), Baskharone and Hensel (1991a, b), and, 

Athavale et al. (1994) use time dependent 3-dimensional computational fluid 
dynamic models to obtain rotordynamic force coefficients for specified small 
amplitude whirl motions about a journal equilibrium position. The intensive 
unsteady flow computations are not coupled to the dynamics of a rotating system, 
and thus, can not yet address effectively the nonlinear aspects of the fluid 
film bearing element studied. Padavala (1993) recently completed a time 
dependent analysis for a rigid rotor supported on turbulent flow annular 
pressure seals. The numerical integration of the unsteady bulk-flow equations is 
performed in explicit form, and thus, the method requires very small time steps 
to assure stability of the calculated solution. Padavala concludes that, for 
most loading cases of interest, annular pressure seals in cryogenic 
turbomachinery behave as linear elements even for large amplitude motions away 
from the rotor equilibrium position. The author also notes that fluid inertia 
forces are unimportant in the unsteady bulk-flow model. However, Padavala fails 
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to indicate that the test cases studied are relevant to a rotor system operating 
well below its critical speed. This operating condition alone justifies his 
major conclusion. 


ANALYSIS 

Consider a point mass rigid rotor supported on fluid film bearings as shown 
schematically in Figure 1. The equations describing the motion of the 
rotor-bearing system are given as: 


M A x = F x (t) + V x (t) + M-g 
M A y = F y (t) + V y (t) 


( 1 ) 


where {A x ,A y } correspond to the rotor accelerations in the X and Y inertial 
coordinates, and ( F l > w l ) x=Xf Y corres P ond to the time dependent fluid film 
bearing reaction forces and applied external loads, respectively. The solution 
of equations (1) is unique provided initial conditions on the rotor center 
displacements {X o ,Y q } and velocities {V Xo ,V yQ } are specified at some initial 
time (t^). In general, the fluid film bearing forces depend on the bearing 
geometry, operating conditions and the journal (rotor) center displacements and 
velocities, and expressed in general form as: 


F x = F i (X,Y,V x ,V y ), i=X, Y (2) 

A local coordinate system is placed on the unwrapped plane of a hybrid 
(hydrostatic/hydrodynamic) fluid film bearing with the {x,y} axes pointing in 
the circumferential and axial directions, respectively (see Figure 2). This 
cartesian coordinate system is a logical consequence of the smallness of the 
film thickness relative to the bearing diameter and length and constitutes the 
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fundamental assumption of lubrication theory. 

The fluid reaction forces are determined from the film pressure generated 
by the flow field within the fluid film bearing, i.e. 


F = 

i 


2nR 


P h x dx dy, i=X, Y 


(3) 


where h x -cos6, and h y .sin8, respectively. The flow field (velocities and 
pressure) within the bearing is described by a set of equations of notion with 
appropriate boundary conditions and accompanied by specification of the 
kinematics of the journal motion and the film thickness. 

At operating conditions, the journal position relative to the bearing 
housing is described with reference to the inertial axes (X.Y.Z.y) by the 
journal center displacements (X(t),Y(t». The film thickness and its time 
derivative in the flow region are given by the following expressions: 


H(6,y,t)= C(y ) + X(t) cos(G) + Y(t) sin(0) 


9=x/R (4) 


3H/9t = V x cos(0) + v y sin (0) 
where C is the bearing radial clearance, in general a function of the axial 
coordinate. Clearance functions of interest range from uniform and tapered 
(convergent, divergent) to discontinuous functions like steps and more complex 
forms of the structural wavy type. The film thickness above is representative of 
a cylindrical bearing. The analysis can also represent film thikness variations 
due to pad preload in fixed arc or tilting pad bearings, and due journal 
misalignment axis conditions. In equation (4), the journal center displacement 

components (X.Y) are functions only of time and determined by solution of the 
governing equations (1). 
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T he equations of motion for the fluid film bearing 

The unsteady fluid flov on the film lands of a bearing is considered as 
fully developed with a turbulent character due to the large axial pressure drop, 
high rotor surface speed and the lov viscosity typical of process liquids. Here, 
"average fluid inertia" equations are used to describe the bulk-flov motion in 
thin film geometries (Elrod et al. , 1983, Sreri et al„ 1983). The equations of 
mass conservation, momentum and energy transport for the turbulent bulk-flov on 
the thin film lands of a bearing are given by Yang et al. (1992 , 1993a, 1994) asi 
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i=x,y 

vhere A^.O-R and A y ,0 correspond to the surface rotor speed in the circumfer - 
ential and axial directions, respectively. V K y .( Kj „ B ), 2 are the vail shear 
stress difference coefficients as local functions of the turbulent friction 
factors, bearing and journal surface conditions, and the flov Reynolds numbers 
relative to the rotating ( Rj) and stationary (Rg) surfaces, i.e V f jRj , k^R, 
(Hirs, 1973). For cryogenic liquids such as LH 2 , LO.,, LN^, and LCH 4 , the fluid 
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properties are calculated from the Benedict-Web— Rubin equation of state as given 
in the standard data base of McCarty(1986) . And, 

V h B< T - T B> * h J (T - T J> (8) 

denotes the heat flow from the fluid film to the bearing and journal surfaces at 
temperatures (Tg) and (Tj), respectively, and with (hgjhj) as the bulk-flow 
convection heat transfer coefficients. A description of these parameters in 
turbulent bulk-flows is given by Yang (1992). 

In an orifice-compensated hydrostatic bearing, mass conservation must be 
satisfied at each bearing recess of area (1-R-0J and depth H r . The global 
balance between the mass flow through the orifice restrictor (Q ro ), the mass 
flow into the film lands and the accumulation of mass within the recess volume 
(V ) is given as: 


Q , =A J 2 p (P -P " ) 
rk o H r v s rk' 


-» -+ 
p H U.n 


3V , r 9P 1 3T . I 

~ Vf ] <»> 


‘ rk 


for k=l,2..., Nrecess 

2 

where A o =C d nd /4 is the effective orifice area, and P p =(l/p) 3p/9P and 
Rj*=-( 1/ p) 3p/ 9T= represent the fluid compressibility and volumetric expansion 
coefficients, respectively. T r represents the closure (perimeter) of the recess 
volume with the film lands and with normal n along the boundary line. Note that 
the orifice flow equation is valid only for small changes of the liquid density. 
A more accurate relationship for highly compressible cryogenic liquids is not 
readily available, although a complex thermophysical model has been advanced by 
Hall et al. (1986). Experimental evidence has shown that the predicted 
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performance of hydrostatic bearings is quite sensitive to the appropriate 

selection or knowledge of the discharge coefficient (c,). See Scharrer et al. 

(1990), Kurtin e, al. (1993), and Franchek et al. (1994) for more details on 
this subject. 

The fluid recess-edge pressure at the entrance to the film lands is given 
by the superposition of viscous shear effects on the recess extent and an 
entrance drop due to fluid inertia. Figure 3 shows the assumed pressure 
distribution within the recess volume and details the relevant nomenclature. On 
the circumferential direction, the pressure rise (P*) downstream of the recess 
orifice is given by (Constantinescu and DiMofte, 1987, San Andres, 1992): 





R-e 

p e = P r - UK ? 
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U x (p e /p r ) 11 - ' 
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2 H Z 
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( 10 ) 


k=l,2. . , Nrecesses 

where, H is the circumferential flow local Mach number at the orifice discharge 
and defined as the ratio between the azimuthal velocity U,, and the liquid sound 
speed (l.//[SpP r I). The entrance pressures (P*) to the film lands in the 
circumferential and axial directions are given by: 


*1 • p e' - r*<i*v 1- V K 


i=x,y 


( 11 ) 


for k*l,2. . . , Nrecess 

The analysis generalizes equations (11) for uneven empirical entrance loss 
factors 5 in the axial direction (y) and also circumferentially(x) upstream(u) 
and downstream(d) of the recess. The Bernoulli like pressure drop occurs only if 
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fluid flows from the recess towards the film lands. On the contrary, if fluid 
enters from the film lands into the bearing recess, then the edge pressure is 
equal to the recess pressure. This consideration is based on momentum 
conservation for turbulent shear flows in sudden expansions. 

A simplified form of the energy transport equation at the recess volume is 
easily determined from equation (7). A detailed exposition of the energy 
transport process in a bearing recess is given by San Andres (1993). 

Boundary Conditions 

At low rotational speeds, the pressure at the bearing side boundaries is 

invariant in time and equal to a specified value of ambient or sump pressure, 
i . e. , 


P (*,0) = P Q ; P (x,L) = P L (12) 

For multiple-pad bearings, at high journal surface speeds (2-R) significant 
momentum changes at a pad leading edge generate a "ram" pressure which is a 
fraction of the dynamic head based on the journal surface speed (Burton and 
Carper, 1967, Ha et al., 1994). Details of the analysis including the inlet 
fluid speed and energy transfer at the leading edge of a bearing pad are given 
in a previous publication (San Andres, 1993). The axial velocity U y is zero if 
the bearing is symmetric both in geometry and operating conditions. The symmetry 
condition allows simplification and solution of the flow field in one-half of 
the bearing. 

Numeric al Method of Solution for the Fluid Flow Equations 

A cell finite-difference scheme is implemented to solve the nonlinear, time 
dependent differential equations on the film lands. The numerical procedure is 
based on the forward marching scheme of Launder and Leschziner (1978) and uses 
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the SIMPLEC algorithm of Van Doormaal and Raithby (1984). At a bearing recess, 
the flow through the supply orifice is coupled to the flow towards the film 
lands and to the unsteady accumulation of fluid volume on the recess. A 
Newton-Raphson scheme is then used to update the recess pressures and to satisfy 
the transient mass continuity constraint at each bearing recess. 

Appendix A details the algebraic difference equations derived on a 
staggered grid for calculation of the flow velocities and pressure on the 
bearing film lands. The time marching algorithm is fully implicit with the flow 
variables (velocities, pressure and temperature) obtained from numerical 
solution of the equations of motion at the current time step and accounting for 
the past history of the flow field. The scheme is computationally intensive and 
requires the satisfaction of boundary conditions and reccess mass flow equations 
at each time step. The precursor numerical procedures to the present time - 
dependent implicit scheme have been extensively validated for steady state flow 
problems with existing measurements of bearing static load performance and 
rotordynamic force coefficients (San Andres, 1990, 1992). Kurtin et al.(1993), 
Francfiek et al.(1994), Mosher (1992) present detailed experimental vs. 
theoretical predictions for water lubricated, turbulent flow HJBs. Further 
validations to experimental force coefficient data for LH 2 HJBs and water HJBs 
are reported by Yang et al. (1994). San Andres (1993) also presents comparisons 
to test force coefficients for an oil lubricated HJB from Adams et al. (1992), 
journal bearing load measurements from Tonnesen and Hansen (1981), and numerical 
predictions by Braun et al. for the performance of liquid oxygen journal 
bearings (1987). 

In general, computations for journal equilibirum position at steady state 
operating conditions require a relatively small number of grid points to achieve 
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grid independent results. For example, less than 2 % difference in bearing static 
and dynamic performance characteristics are determined when comparing the 
results from a 67 by 16 grid (number of circumferential points x axial points) 
with those from a 85 by 26 grid for bearing cases reported by Yang (1992). Note 
that the analysis and program developed can handle also fluid film bearings for 

industrial applications and operating in the laminar-flow regime with viscous 
lubricants . 

Several empirical parameters are used in the analysis of hydrostatic 
bearings and annular pressure seals. Entrance loss coefficients (^,^) and 
rotor( journal) and stator(bearing) surface roughnesses (r^r^ are needed for 
annular seals and HJBs, while orifice discharge coefficients (C^) are also 
required for hydrostatic bearings. Lindsey (1993) presents a detailed 
sensitivity analysis for annular seals, and Kurtin et al. (1993), Mosher (1992), 
and Franchek et al. (1994) present similar studies for hydrostatic bearings. 

These references describe comparisons between numerical predictions based on the 
present model and experimental measurements performed on water lubricated seals 
and hydrostatic bearings at an existing test facility (Childs and Hale, 1994). 

The sensitivity analysis involved changing an input parameter by ±10% from its 
estimated experimental value for each operating condition while the other input 
parameters were kept invariant. The maximum difference between the numerical 
prediction and experimental value for each case was then compared with the 
maximum difference from the original results. Then, a relative sensitivity was 
determined by dividing the percentage change in maximum difference of the 
calculated parameter by the percentage change in the input parameter. The 
studies determined that flow rate and direct stiffness coefficients are 
particularly sensitive to changes in the orifice discharge coefficient (C^) for 
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HJBs and less sensitive to variations in the entrance loss coefficients. On the 
other hand, flow rate and direct stiffness coefficients for annular pressure 
seals are highly sensitive to the changes in the inlet loss coefficient (^) and 
surface roughness conditions. Other force coefficients and the whirl frequency 
ratio for HJBs are particularly insensitive to variations in the empirical 
parameters. The interested reader is encouraged to seek the references cited 
for further details on the sensitivity analyses and major results. 

Nu merical Solution of the Equations of Hotion for the Coupled 
Rotor-Bearing System 

The general nonlinear equations (1) governing the rotor and bearing motion 
are numerically solved using the Vilson-0 method as described by Bathe (1982). 
Selection of this implicit numerical method is based on its stability 
characteristics and accuracy, in particular for linear systems. A brief 
description of the method as applied to the present analysis follows. 

The time domain is divided into a sequence of discrete times ft 1 

j j-1* Nt 

starting at the initial time t Q , and such that the time step Atj is related to 
the discrete times by (t. +1 -t.). At discrete times t. and t.=t. +r the equations 
of motion (1) are reduced to the following algebraic expressions: 


M A Xi ~ F Xi + W Xi + M ' g 
M A Yi = F Yi + W Yi 


MA Xj " F Xj + W Xj + M ’ g 
MA Yj = F Yj + W Yj 


(13) 


where, A^^ ^(t^), anc * ^Xi - ^X^i^ are t ^ le ^-direction rotor acceleration and 
external load evaluated at time t.. F xi = F X (X.,Y., V x ., Vy.) is the X-fluid 
film bearing force estimated at time t.. Similar definitions hold for the Y- 
direction journal displacement and its time derivatives, as well as forces at 
the calculation times t. and t . . Note that in the equations above, the journal 
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position and velocity and the fluid film forces are to be determined at time t., 

J 

and henceforth all quantities at time t. are known. Equations for incremental 
changes in the system acceleration are determined from equations (13) as: 


M M Xi = AF Xi + * u xi 
M M Yi = AF Yi + AW Yi 


(14) 


where AAX.= Ajy-A^, AVX.= V^., AF xi = F xj -F xi , etc. are referred as the 
incremental changes in accelerations, loads, and bearing forces, respectively. 

For small time steps At^ the fluid film forces using a Taylor series expansion 
are approximately given as: 
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where the local or linearized stiffness, damping and inertia coefficients for 
the fluid film bearing are given as 


3F 
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0 1 (Xi , Yi ) 
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(X. ,Y. ) 

l l 


M a0 , ; «,e=x,Y (16) 

3A B IWi.Tj) 


and evaluated from analytical perturbations to the nonlinear flow governing 
equations for the fluid film bearing. Complete details on the derivation and 
solution for the flow equations for small perturbations about an operating 
condition have been given in detail in the past (San Andres, 1993, Yang, 1992). 
In general, for turbulent-flow bearings and seals with compressible fluids, the 
linearized force coefficients are functions of the excitation whirl frequency. 

Substitution of the force coefficients into equations (14) for the 
incremental accelerations leads to the following linearized equations: 
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'* M xx m xy jr M xi| + Fxx K xy]Fxi] ^ T c xx c xy] ['“i] p w K i 

Hyx M*M Yr J[M y .J Ik yx K yy J U y J [cjjj C yy J U. J " L„, 


In the Wilson-0 method the accelerations are regarded as linear on the 
extended time interval t. to t. + 0At.. The method depends on the parameter 0 with 
optimum results for a value 0=1.42 (Bathe, 1982). Thus, over the local time 
coordinate t the accelerations are given as (see Figure 4): 

A x (t) = AX. + (x/S.) AsA x . 

A y (T) = A + (t/S.) AsA y . , S.«0At. Vl8> 


1 - 1 


Integration of these equations over the discrete time 0<x<S i gives the 
displacements and velocities as quadratic and cubic functions of time x as: 

V X (T) = V Xi + T A Xi + ( t2/2S i) A s A xi 

and, " y(T) = + T "Yi + (t2/2 V Vli (19) 

X(T) = X. + V x . X + (t 2 /2) A xi + (x 3 /6S.) AsA x . 

Y(x) Y i + V Y i X + ( T /2) A yi + (t 3 /6S.) AsA y . (20) 

Substitution of x=S. on the expressions above leads to relationships for the 

incremental changes in acceleration and velocities at the end of the extended 
time interval, i.e. 


AsA Yi 


(6/S.) AsX. - (6/S.) V x . - 3 A x . 
(6/S 2 ) AsY. - (6/S.) V y . - 3 A 


AsV 

Xi 

AsV 

Yi 


(3/S.) AsX. - 3 V xi - (Sj/2) A x . 
(3/S.) As y . - 3 V y . - (S./2) A y . 
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From equations (17), the equations for evaluation of the incremental 
accelerations, velocities and displacements at the end of the extended time 
interval are given as: 


* +M xx m xy 
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(23) 

The expressions for the incremental changes in accelerations (AsA x y ) and 
elocities (AsV^ y ) are substituted into these equations to give: 


[Hi) {(6/S?) AsZ. - (6/S.) V zi - 3 } 
+ [Ci] {(3/S.) AsZ. - 3 V Z1 - (S./2) A zi ) 
+ [Ki] AsZ. = 0 { W(t j )-W(t.) } = AsV. 


where [Ki], [Ci], [Mi] represent the matrices of local stiffness, damping and 
inertia coefficients evaluated at time t.. AsZ = {AsX, AsY} T is a vector of 
incremental displacements, and AsV = [0-AV^ Q.AV Y J T is a vector of incremental 

external loads. Rearranging equations (24) gives finally the following matrix 
expression: 



[K ] AsZ. = AsV? 

i l 

(25) 

where, 

[K*J = [Ki] + (3/S^ ) [Ci] + (6/S?) [Mi] 

(26) 

and 

AsV* = ©jv.-V.j + (6/S.)[HiJ + 3 [ Ci ] Vj,. + 

3[Mi]+(S./2)[Ci] A. 


are known as the "pseudo" dynamic stiffness matrix and equivalent load vector 
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respectively. Note that equation (25) represents the dynamic problem as a 
"pseudo" static load problem (Bathe, 1982). 

Once the incremental displacements over the extended time period S. are 

obtained from solution of equation (25), the following formulae are used to 

estimate the displacements and velocities at the end of the time step t.=t + At 

J i i 

V Zj = V Zi + Ati t A zi + ( 1/2 > ** zi ) 

Z j = Z 1 + At i V z . + (AtJ/2) { A z . + (1/3) AA^ } <28> 

where r 

M z . = ( 1 / 0 ) AsA z . = 1 / 0 J ( 6 /S?) AsZ. - ( 6 /S.) V z . - 3 A z± j (29) 

To insure dynamic equilibrium, the accelerations at the end of the time step 

(tj > are not determined from the incremental equations but by direct solution of 
the original equations (1), i. e . 

A Zj = < 1/M > t F Z j + + M ’*J (30) 

where the fluid film forces are calculated from the solution of the flow 
equations with the updated values of journal displacements and speeds on the 
film thickness and its time derivative, i.e., F = F (Z V ) 

At the start time (to) for analysis, initial conditions in journal position 
and speed and the corresponding flow field (pressure, temperature, velocities 
and fluid material properties) must be given. Then, for time specified 
variations in the external loads, the procedure follows the scheme outlined 
above at each time step. A User's Manual describes the operation of the computer 

program developed along with the many options available and bearing geometries 
analyzed (San Andres, 1994a). 
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Required time step for the transient analysis. 

The computational effort required to solve the entire bearing flow field along 
with evaluation of forces and local force coefficients is very intensive. It is 
noted that time steps need to be kept small if accurate results are to be the 
natural outcome of the numerical algorithm. The minimum required size of the 
time step is typically related to the fundamental period of motion of the system 
or the largest frequency intrinsic to the external load. Thus, before 
calculations are started an estimate of the natural frequency of the system for 
small amplitude motions about an equilibrium position is required. Note that the 
notion of a natural frequency here implies linearity of the equations of motion 
which is not likely to be the same as in a highly non - linear system with large 
time varying forces applied to it. 

The criteria for the threshold speed of instability (Q s ) and whirl 
frequency ratio (♦) of a rigid rotor supported in fluid film bearings are well 
known (Lund, 1965). For fluid film bearings where fluid inertia effects are not 
important the equivalent stiffness and whirl frequency ratio of the rigid rotor 
- bearing system are given as: 



K XX C YY +C XX K YY~ C YX K XY~ C XY K Yy 


C XX +C YY 


(31) 


* 2 = (« n /sy 2 = (Ke q~^ xx)(Ke q" KYY) ~ KxYK Yx 

Q s ^ C XX C YY -C XY C YX ^ 


(32) 


where the force coefficients {K^.C^} are evaluated at an equilibrium journal 
(rotor) position defined by the applied static load. From these equations, an 
estimate of the natural period of motion (xj of the rotor-bearing system is 
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given as: 


T = 2 n/w 
n n 


2n (M/K ) 
eq / 


1/2 


(33) 


If the operating speed (0) is larger than Q s then the rotor bearing system 
becomes unstable (in the sense of a linear stability analysis) vith unbounded 
motions of the rotor at its natural frequency «»„), i.e. at a fraction * of the 

rotational speed. For a similar analysis including fluid inertia effects refer 
to San Andres (1991). 

As briefly stated earlier, the time step for the transient analysis must be 
a fraction of the natural period of motion (T,), or a fraction of the period of 
motion due to rotor angular speed (2n/S), „ r a fraction of the smallest period 
determined by the largest frequency («,) intrinsic to the time-varying external 
loads. The time step selected should be the smallest satisfying these three 

conditions. In practice, the time step is chosen for purposes of accuracy as 
(Childs et al. , 1977) : 


At < [2it/w ]/60 


(34) 


there is the largest frequency (in Hz) likely to be excited in the system. 
The less stringent requirement, At < 2/ % , guarantees only stability but not 
necessarily accuracy of the numerical scheme. Note that suddenly applied loads 
like impulses and step loads are rich in high frequencies, and consequently 
initial time steps are constrained to be very small in order to capture vith 
accuracy the imposed external actions. 

The efficiency of the complex computational model coupling the dynamics of 
the rotor-bearing system to the flow dynamics within the fluid film bearing is 
severely reduced due to the stringent requirement of small time steps and the 
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current limitations in the processing speed of computers. The current program 
runs on a personal workstation with a clock speed of 35MHz. Times of execution 

reported later for the examples and test cases analized are particular to the 
machine used. 

For example, calculation of the steady-stated flow field and forces in a 
water lubricated hydrostatic bearing operating at 30 kcpm and given specified 
static load is about 60 seconds with an initially assumed flow field close 
within 20 X to the actual solution. This execution time seems small (in my 
opinion) if a few static load cases are analyzed. In the past this fast 
execution has allowed the current computational program to be used in an 
interactive way with the user. However, to describe a transient situation due to 
a blade loss (sudden unbalance) turns out to be almost prohibitive. If we assume 
that the entire transient will take about 10 rotor revolutions then the total 
time of calculation is equal to 10 x 0.002 secs where this last number 
corresponds to the period of whirl motion (1/500 Hz). With a time step 
At=0. 002/60 (33.3 usees) to complete the necessary 600 time steps, the present 
computational model would approximately need about 600x60=36,000 secs (10 hours) 
of computing time. This time value is certainly excessive and addresses to the 
need of either refining the analysis, coming across with more efficient ways to 
calculate the system transient response, or simply getting hold of a faster 
computer. 

The recent archival literature presents several transient analyses for 
rotors supported on inertialess - fluid film bearings, see for example Padavala 
(1993), Gadangi et al. (1994a, b), Desbordes et al. (1994). These references 
report advanced models for the prediction of fluid film bearing forces, but 
however little insight is given in regard to the time steps used and the actual 
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computational time required to complete a detailed transient analysis. The 
authors of these scholarly references assert in personal communications that 
their schemes are time consuming and need to run on the fastest computers 
available, namely a CRAY machine, for a few hours. 


Approxi mate (linear) transient solution 

Externally pressurized bearings like hydrostatic bearings and annular 
pressure seals are know, to provide force coefficients which do not vary greatly 
with the journal position (or the applied load) at a fixed operational speed 
(San Andres, 1991). This is not the case, however, for hydrodynamic journal 
bearings as the technical literature profusely reports. The uniformity of force 
coefficients in HJBs and pressure seals allows a dramatic simplification of the 
general nonlinear equations of motions (1) to render: 


(35) 


° r lMoJ M Z + t C °] AV Z + IK°] ^ = tf( t ) (36) 

AX-X-Xo, AV X =V X , M x =A x , etc. represent the approximate journal (rotor) 
center displacements, velocities and accelerations about the (equilibrium) 
steady-state position (Xo,Yo) defined by the static load (M-g). The set of 
f lcients {K^, C^, M^} a> p =x>Y correspond to the linearized stiffness, 
damping and fluid inertia force coefficients evaluated at (Xo,Yo). The numerical 
solution of the equations above is much faster than the full nonlinear solution 
described earlier since the matrices of inertia, damping and stiffness are 
evaluated only once prior to the start of the transient analysis, and 
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furthermore, the method does not even require the solution of the flov equations 
at each step. The reader is urged to note that the approximate linear equations 
of motion become exact for small amplitude journal motions induced by small 
magnitude loads, and in effect, under these conditions they can even be used for 
rotors supported on hydrodynamic journal bearings. In fact, the greatest 

achievements of rotordynamic analysis have become fro. such a simple and -crude" 
representation of the system dynamics. 

The method selected to solve the approximate linear equations is the 
Wilson-0 scheme described earlier. Execution times for solution of the 
approximate linearized equations of motion are about 4 orders of magnitude lower 
than the full nonlinear transient solution. Examples detailed below compare the 
relative accuracy of the approximate analysis when correlated directly to the 
full nonlinear transient model. 

RESULTS AND DISCUSSION 

The computer program named hydrotran is the outcome of the fluid film 
bearing transient analysis. Details on the capabilities and major options 
available with this program are fully described on its Manual (San Andres, 

1994a). All the transient cases studied including the input and output files as 

well as the graphical output are given on the program Examples (San Andres, 

1994b). 

This section discusses the transient response of several point mass rotors 
supported on different types of fluid film bearings like plain cylindrical 
bearings, annular pressure seals and hydrostatic bearings, and under a variety 
of impossed external loads such as step forces, Impulsive loads and suddenly 
applied unbalance forces. First, the examples compare the calculated transient 
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responses with those given in the literature for some relevant cases. Next, 
results are presented for a cryogenic application with annular pressure seals 
(damper seals) and also the time response of a high speed, water lubricated 
hydrostatic bearing and rotor system. No efforts have been made to compare the 
numerical predictions to experimental results since there are none reproted in 
the archival literature. 

Table 1 presents the description of a rotor supported on plain journal 
bearings as given by Childs et al. (1977). The example presents the transient 
response of the system due to a step load equal to the rotor weight (8 kN). The 
journal (shaft) operates at a speed of 6,493 rpm and it is released initially 
from a position at the center of the bearing clearance. Figure 5 shows the 
the journal center to move rapidly from the centered position towards the bottom 
of the bearing clearance. At a large eccentricity the journal suddenly reverses 
its motion and proceeds steadily to its equilibrium static position. The results 
depicted in the figure are very similar to those of Childs et al. (1977) and 
also confirmed with a prior transient response on short journal bearings (San 
Andres, 1984). Figure 6 shows the transient fluid film forces (F V ,F V ), rotor 
inertia forces (M*A x ,M*Ay) due to the external loads (W x =8 kN, Vy= 0 N). The 
numerical results feature very large inertial rotor forces and bearing reaction 
forces in the direction of the aplied load (X) at the instant of journal 
reversal motion. The maximum bearing X-force is actually 7 times larger than the 
applied static load. 

Childs et al. (1977) present an interesting example for a rotor bearing 
system with a highly nonlinear behavior. Table 2 describes the geometry of the 
rotor and hydrodynamic journal bearings. At the operating conditions the journal 
is initially located at its equilibrium position due to the rotor weight. Then, 
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a sudden unbalance load (pressumably due to a blade loss) of a magnitude equal 
to 697 N is exerted on the system. Numerical predictions for the bearing rotor- 
dynamic force coefficients at the static equilibrium condition show that the 
rotor-bearing system is linearly unstable at the rated speed. Figure 7 shows the 
journal orbital motion to proceed from the static position to describe a large 
orbit within the bearing clearance. The journal motion is very similar as that 
reported by Childs et al. (1977). Figures 8 and 9 show the "steady state- 
journal displacements and bearing reaction forces after 10 rotor revolutions 
have elapsed, i.e. for times greater than 57.1 msec. The "T" labels on the 
figures mark a full revolution of the rotor. The numerical predictions show a 
characterisitic half frequency whirl motion and demonstrate a highly nonlinear 
behavior. Note in particular the complex functional form of the bearing reaction 

forces which have a period twice as large as that of the impossed unbalance 
force. 

Tichy and Bou-Said (1991) present a study on the effect of impulsive loads 
on rotors supported on -short length" plain journal bearings. The thrust of the 
analysis is to demonstrate the importance of fluid inertia effects on the 
initial time response of rotor-bearing systems. Table 3 presents the geometry 
and operating conditions of the rotating system studied by Tichy et al. At the 
static equilibrium position (due to the rotor weight) the system is linearly 
stable. Note that the current analysis predicts a different journal equilibrium 
position if fluid inertia effects are included in the flow model equations. The 
impulsive load has an exponential character and decreases to 36.7* its initial 
value after t Q =2.77 msec have elapsed (1/6 rotor period). At the time the rotor 
has given a full revolution (t=16.6 msec), the impulse load has practically 
dissapeared. Figure 10 shows the journal orbital motion for the applied 
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impulsive load along the -Y axis, and Figure 11 depicts the journal center 
isplacements (X,Y) for the time span of numerical integration. Figures 12 and 
13 show the bearing reaction forces and rotor inertia forces for the initial 
transient motion to 40 msec (approximately 2.40 rotor revolutions) , respectively . 
The inertialess flow model "appears" to predict a more rapid decay response than 
the fluid inertia model; however, note that both flow bearing models approach 
different equilibrium positions. The inertial flow model does not predict the 
"ripple" at the end of the journal trajectory in the direction of the applied 
impulsive load. This condition may be solely due to the smallness of the time 
step used. The most important conclusion from the example is derived from the 
i nitial bearing reaction forces and rotor mass inertia forces which are 
substantially smaller for the fluid inertia model. The model without fluid 
inertia predicts bearing forces up to 66* larger than those from the fluid 
inertia model. The reason for the difference lies on the "added" bearing inertia 
forces which arise due to fluid inertia effects. It is relatively simple to 
demonstrate that the apparent mass of the fluid is about 40 kg, and determines 
an overall increase in the system inertia to 245 kg. Hence, the predicted 
journal accelerations are smaller for the fluid inertia model, and consequently 
the rotor inertia force will be substantially reduced. The present transient 
response results are somewhat different to those reproted by Tichy and Bou-Said 
(1991). The reference presents an orbital plot where the fluid inertia solution 
shows the "ripple" at the farthest position of the journal trajectory. 

Furthermore, even their steady-state equilibrium results could not be verified 

with the exact formulae for the inertialess finite-length bearing model given by 
Lund (1966). 

In the examples for hydrodynamic journal bearings above no comparisons with 
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the approximate (linear) solutions are presented since these types of fluid film 
bearings are well known to be highly nonlinear. Furthermore, the magnitude of 
the external loads applied are large enough to prevent the use of the linear 
transient solution. The Examples Handbook (San Andres, 1994b) does include a 
comparison of linear and nonlinear rotor-plain bearing transient responses using 
a test case from Gadangi et al. (1994a). 

Table 4 presents a description of a cryogenic rotor-bearing system. The 
fluid film bearings are tapered damper seals handling liquid oxygen. Predictions 
for the static and rotordynamic force coefficients of this type of seal are 
given in detail by San Andres and Yang (1994). This example compares the 
transient system response obtained using the nonlinear and linear time 
integration models. The seal rotordynamic force coefficients at its centered 
position are used in the linear model. The predicted "linearized" natural 
frequency of the system is 290 Hz, and therefore, the rotor speed of 26,000 rpm 
is 1.5 times greater than the system critical speed. 

Figures 14 to 16 present the seal nonlinear and linear transient responses 
for a step and a ramp loads applied in the X-direction. The magnitude of the 
load is 10 kN with a rise time equal to 11.5 msec for the ramp load. The orbital 
seal motions are shown in Figure 14, and the transient seal eccentricy is 
depicted in Figure 15. The results show a forced motion with little damping and 
at a frequency close to the predicted natural frequency for small amplitude 
motions. Note that the linear (approximate) transient solution follows the same 
trend as the nonlinear solution but presents a smaller overshoot and decays 
faster. The transient bearing forces (F x ) in Figure 16 show also a system with 
very little damping. For the step load the initial overshoot is large 
(approximately 100pm or 77% of the minimum clearance). Even for this large 


120 



amplitude motion, the approximate solution predicts well the system response 

since the damper seal has force coefficients which vary very little with the 
journal position. 

Figure 17 shows the rotor-damper seal transient response to an impulse load 
of magnitude equal to 10 kN and with a characteristic time of 0.57 msec. The 
seal reaction force (F y ) in the direction of the applied impulsive load (W y ) is 
depicted in Figure 18. Note that after a full rotor revolution (2.3 msec), the 
applied impulse vanishes. The linear (approximate) solution predicts well the 
initial seal forces but it presents too much damping and decays faster than the 
nonlinear transient response. 

The time response of the rotor-damper seal to a synchronous periodic load 
is presented next. The magnitude of the load is 5 kN which corresponds to a 
rotor unbalance of 13.5 pm. Figure 19 shows the transient seal eccentricity and 
Figure 20 depicts the magnitude of the damper seal reaction force for the linear 
and nonlinear models. Figure 21 shows the transient seal center orbital motion 
for the nonlinear model. The linear model predicts well the transient response 
but again it shows more damping and reaches faster a steady-state solution which 
is slightly larger than that of the nonlinear model. The correlation between the 
linear and nonlinear transient models for the example rotor supported on liquid 
oxygen damper seals is regarded as satisfactory and demonstrates that externally 
pressurized bearings can be well represented by linearized models. 

Table 5 presents a description of a rotor supported on water lubricated 
hydrostatic bearings (HJBs). The bearings have five rectangular recesseses and 
orifice ports evenly distributed around the bearing circumference. Kurtin et al. 
(1993), Childs and Hale (1994), and Franchek et al. (1994) report experimental 
results for the bearing flow rate, load capacity and force coefficients of this 


121 



bearing tor verious supply pressures and speeds. The tests show that the bearing 
is very stiff at high supply pressures and demonstrate that the force 
coefficients have a weak dependence on the journal position. Linear and 
nonlinear transient responses are compared in the discussion that follows. The 
bearing force coefficients at the journal centered position are used in the 
approximate (linear) transient model. 

Figures 22 and 23 show the hydrostatic bearing transient response to 
applied (5 kN) ramp loads with rise times equal to 0.2 msec and 2 msec, 
respectively. Note that the rotor period of motion is 2 msec, and the fast-ramp 
load reaches its maximum value in 1/10 of a rotor full revolution. Figure 24 
depicts the bearing reaction force (F x ) in the direction of the applied ramp 
load. The numerical predictions show the system to be highly damped with similar 
responses for both the linear and nonlinear transient models. 

Figures 25 and 26 show the transient response of the rotor-HJB system to an 
impulse load of 5 kN with a characteristic time of 0.5 msec. The bearing 
reaction force (F Y ) in the direction of the applied impulsive load (W Y ) is 
depicted in Figure 27. Note that after a full rotor revolution (2 msec), the 
applied impulse vanishes. The linear (approximate) solution shows slightly more 
damping than the full nonlinear model. Both models compare well in all aspects 
of the time response since the maximum journal displacements are small relative 
to the bearing clearance (76.2 pm). For example, Figure 26 shows that the 
maximum journal center displacements are 9.8 pm and 9.5 pm for the nonlinear and 

linear models, respectively. These values represent motions to just 13* of the 
bearing nominal clearance. 

The time response of the rotor supported on hydrostatic bearings due to 
rotor unbalances is presented next. Two unbalance displacements of 15.2 pm and 
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40 pm determine loads of magnitude equal to 7,520 N and 19,740 kN, respectively. 
Figure 28 shows the transient journal eccentricity for 15 rotor revolutions, and 
Figure 29 depicts the steady-state rotor orbits. For the largest unbalance (U=40 
um), the maximum journal displacements reach to 90% of the bearing clearance. 

The steady-state orbit magnitudes for the linear model are different (within AX 
of the bearing clearance) to those predicted by the nonlinear model. For the 
largest unbalance, initial and steady-state bearing reaction forces and rotor 
inertia forces are presented in Figures 30 and 31, respectively. The linear and 
nonlinear models correlate well and show that operation at 30 kcpm is above the 

system natural frequency since the rotor inertia forces are almost in phase with 
the unbalance load* 

The linear (approximate) transient model appears adequate to predict the 
transient response of point mass rotors supported on externally pressurized 
bearings (damper seals and hydrostatic bearings). However, the most important 
conclusion is related to the dramatic differences in computing time for the 
linear and nonlinear transient models. For example, the nonlinear (large) 
unbalance response of the example rotor-hydrostatic bearings takes 38,380 
seconds to compute 450 time steps (85.3 sec/step) in a SGI 35 MHz worksattion, 

while the approximate solution takes only about a second. The time difference is 
over four orders of magnitude. 


CONCLUSIONS 

A detailed analysis for the prediction of the transient response of a rigid 
rotor supported in fluid film bearings is presented. The nonlinear equations of 
motion of the rotor-bearing system are numerically solved using the Wilson-0 
method with local space linearization at each time step. The fluid film bearing 
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forces are evaluated from the numerical solution of the unsteady bulk-flow 
equations on the film lands and with satisfaction of the time-dependent flow 
constraints at the recesses of a hydrostatic bearing. An implicit computational 
fluid mechanics algorithm is implemented to solve the bearing bulk-flow 
equations. Analytical perturbations renders linear bulk-flow equations for 
estimation of linearized force coefficients. The current numerical model is 
largely based on efficient and accurate models deviced in the past for the 
prediction of the static and rotordynamic force response characteristics of 
turbulent flow fluid film bearings in aerospace applications. 

Examples for the transient response of rigid rotors supported on 
hydrodynamic bearings, annular pressure seals and hydrostatic bearings under a 
variety of external loads are discussed in detail. Transient responses based on 
an approximate model which uses constant rotordynamic force coefficients are 
also presented. The full nonlinear model demands small time steps and is 

computationally very intensive since it requires solution of the bearing flow 
equations at each time step. 

Externally pressurized bearings generally present rotordynamic force 
coefficients which depend little on the journal equilibrium position at a 
constant operating speed. Hence, the linear transient model should be adequate 
to predict the rotor-bearing system time response due to applied dynamic loads. 
The numerous comparisons between the linear and nonlinear transient responses 
support this assertion. Furthermore, the approximate transient model brings 
dramatic savings on the computing time required to determine the rotor-bearing 
system time response. On the other hand, hydrodynamic bearings are well known to 
be highly nonlinear and the linearized transient response would be accurate only 

for small amplitude loads and at operating conditions well below the stability 
margin of the rotor-bearing system. 
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Table 1. Description of rotor supported on plain cylindrical bearings. 

Transient response due to applied static weight. 

Data from Childs et al. (1977). 

Bearing type: cylindrical plain journal bearing 
C=127 pm, D=50.8 mm, L=25.4 mm, L/D=0.50 
Lubricant: viscosity p=0. 06897 Pa-s, p=890 kg/m 3 
Supply pressure: 0 MPa (0 psig) 

Bearing model: isothermal, incompressible, isoviscous fluid 

Rotor: mass M=816.5 kg, weight M-g= 8,0000 N 

rotational speed Q = 680 rad/s (6,493 rpm). 
period of rotational motion: 9.2 msec 

Initial conditions: X q =Y o =0 m, ^ Q =Vy o =0 m/s [centered position] 

parameters for time integration: At=0. 2msec, Tmax=82.8 msec(8 rotor revolutions) 

46 time steps per rotor period. 

Static equilibrium position for applied weight: X =98.2 pm, Y =53 pm. 

S S 

- o - 

Table 2. Description of rotor supported on plain cylindrical bearings. 

Transient response due to sudden unbalance load. 

Data from Childs et al. (1977). 

Bearing type: cylindrical plain journal bearing 
C=127 pm, D=50.8 mm, L=25.4 mm, L/D=0.50 
Lubricant: viscosity p=0. 06897 Pa-s, p=890 kg/m 3 
Supply pressure: 0 MPa (0 psig) 

Bearing model: isothermal, incompressible, isoviscous fluid 

Rotor: mass M=22.68 kg, weight M-g= 222.3 N 

rotational speed 2 = 1,100 rad/s (10,504 rpm). 
period of rotational motion: 5.71 msec 

o 

Unbalance U: 25.4 pm, force magnitude=M-U- Q = 697 N 

Initial conditions: X q = 3.6 pm, Y q = 19.2 pm, V^ o =Vy o =0 m/,s l e Q u ilibrium position] 

parameters for time integration: At=0. lmsec, Tmax=57.1msec(10 rotor revolutions) 

57 time steps per rotor period. 

Execution time (SGI 35MHz): 2,155 sec (3.8 sec per time step) 

- o - 
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’• fiss-jSi" wrsisjiiir"- 

Data from Tichy and Bou-Said ( 1991 ). 

Searing type: cylindrical plain journal bearing 
C-300 um, D=200 mm, L= 100 mm, L/D=0.50 S 
Lubricant: viscosity u=0.008 Pa-s, p=800 kg/m 3 
Supply pressure: 0 MPa (0 psig) g 

MOdel! !i!; be T 1 : incompressible, isoviscous tluid 
With and without fluid inertia effects 

Rotor: mass M=204.1 kg, weight M-g= 2,007 N 

rotationai speed Q = 277 rad/s (3,600 rpm). 
period of rotational motion: 16.7 msec 

External loads: V x = M -g (rotor weight) 

W Y = - 20 * ex P(-t/t o ) kN, t o =2.77 msec (impulse) 

2 U um t0 v r °T W n ight 

without fluid inertia X°=80.7 pm, pm,’ vfZv^Io m/s 

w AO Yo 

parameters for time integration: At-0 .2msec, Tmax=166msec(10 rotor revolutions) 

oj time steps per rotor period. 

Circumferential fow Reynolds number Re = P 2R C /u -1 111 
Squeeze film Reynolds number - c - U ~ lfl31 > 


Re s =Re c * C/R= 3 ’ 4 


A. equilibrium position: Natural frequency v 35.6 Hz (2,136 rpm) 

Critical mass M =309.7 kg 
system stable c 


o - 
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TabIe 4 ' Sfl'JS.'i"? ° f r ,0 I r Pr, ° rted <l “P« seal bearings. 
Seal data from San Andres and Tang (1994). 

Bearing type: tapered damper seals 

st ( atnr 2 Jh 3 ^ ^>=129.1 D =85.1 mm, L- 22.2 mm, 

stator (bearing) relative roughness = 0.044 (knurled surface) 
Entrance loss factor 5-0.20, Inlet swirl =0.5-9 R } 


Lubricant: liquid oxygen at Ts=110.5°K 
Supply pressure P S : 39.6 MPa (5,750 psia) 
Discharge pressure Pa: 2.09 MPa ( 303 psia) 

Seal model: isothermal, compressible fluid 
with fluid inertia effects 


Rotor: mass M=50 kg, weight M-g= 490 N 

rotational speed fi = 2,723 rad/s (26,000 rpm). 
period of rotational motion: 2.3 msec 


Load cases : 

a) Ramp load : 

b) Step load 

c) Impulsive load: 

d) Periodic load : 


Vx=10 kN-(t/T*), T*=ll . 5 mS ec 
Vx=10 kN for t>0 


Wx=M-g_490 N, Wy=10 kN-exp(-t/to), 
Vx=5 kN-cos(Qt), Wy= 5kN-sin(Qt), 


to=0.57 msec 
Unbalance=13.5 ym 


Initial conditions; 


equilibrium position due to rotor weight 
X q = 2.45 ym, Y q = 0.53 ym, V Xo =V yo =0 m/s 


parameters for time integration: At=33.3ysec, Tmax=46msec(20 rotor 

69 time steps per rotor period. 


revolutions) 


At equilibrium position: 


Natural frequency w = 290 Hz (17,410 rpm) 
system stable n 


Force coefficients at seal centered 

k xx =k yy k xy=- k yx c xx =c yy c xy = ~ c yx 

E-6N/m E-6N/ m E-3Ns/m E-3Ns/m 
17070 34.50 26.63 0794 


position: 

m xx =m yy m xy =_m tx 


Jsg 

0.65 


JS£_ 

0.00 


- o - 
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Table 5. 


Description of rotor supported on hydrostatic 
Bearing data from Childs and Hale (1994). 


bearings. 


Bearing type: orifice compensated, 5 recess hydrostatic bearing 
C=76. 2 pm, D=76. 2 mm, L= 76.2 mm, S 

R ecess dimensions: square 27 mm x 27 mm x 381pm (depth) 
orifice diameter d=1.49 mm, C,=1.0 F ' 

d 


Entrance loss factors ^,^0.0, Inlet swirl =0.5-Q-R 

Lubricant: water at Ts=330°K 

Supply pressure Ps: 6.55 MPa (950 psig) 

Discharge pressure Pa: 0.0 MPa ( 0 psig) 


Model: isothermal, compressible fluid 
with fluid inertia effects 


Rotor: mass M=50 kg, weight M*g= 490 N 

rotational speed 9 = 3,141 rad/s (30,000 rpm). 
period of rotational motion: 2.0 msec 


Load cases : 

a) Ramp load : 

b) Impulsive load: 

c) Periodic load : 


Vx=5 kN-(t/T*), T*=2 msec and 0.2 msec 
Vx=M-g=490N, Vy=5 kN*exp(-t/to) , to*0. 
Wx=Wo • cos( 9t ) , Wy= Vo-sin(St). 

. Unbalance U=15.2 pm, Wo= 7,520 N = 
40.0 pm, Wo=19, 740 N 


5 msec 
M-U-S? 2 


Initial conditions: equilibrium position due to rotor weight 

X q = 0.91 pm, Y q = 1.16 pm, V Xo =V Y q =0 m/s 


parameters for time integration: At=66.6psec, Tmax=30msec(15 rotor 

30 time steps per rotor period. 

At equilibrium position: Natural frequency « n = 264 Hz (15,830 rpm) 

Force coefficients at bearing centered position: 

■WSt K XI- K n C XX’ C YT C xi" -C TJ "xX-'St "xr--"!! 


E-6N/m E-6N/m E-3Ns/ m E-3Ns/m 

211 - 6 315.5 206.5 57.26 


- k 8 
16.66 


- k g 

-2.847 


revolutions) 
stable system 


o 




Figure 1. Rigid rotor supported on fluid film bearings 
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2nR 


Figure 2. Geometry of a Hydrostatic Bearing 
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Conceptual description of pressure rise 
and drop at recess edge of a hydrostatic 
pad bearing, and pressure ram effect at 
leading edge of bearing pad. 
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Acceleration 



Figure 4. Linear acceleration in numerical integration method 






Example from Childs et al. (1977), M=816.5 kg 
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Figure 6. Fluid film forces and inertia forces for transient response 
Sudden step load Wx=Mg=8,000 [N] Wy=0 




Na 
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Example from Childs et al. (1977), M=22.68 kg, 10504rpm, c=127um 



[N] seojoj s-8 Buueeg 
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Figure 9. S-S fluid film forces for transient response. 

Sudden unbalance load Wx=Mg+M u co2 cos(cot), 

Wy=M u co2 sin(cot); u=25.4um 




Example from Tichy & Bou-Said- M=2 
0.0003 i t 



Na 
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X [m] 

Figure 10 Effect of fluid inertia on 
Impulsive load Wx=M.g=2KN, W' 
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time (secs) 

Figure 12. Bearing forces with and without fluid inertia 
Impulsive load Wx=M.g=2KN, Wy=-20KN over 2.77msec; 



Example from Tichy & Bou-Said: M=204.8kg,3600rpm,C=300um 



[jsj] seojoj eipeu! jotoy 
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time (secs) 

Figure 13. Rotor Inertia forces with and without fluid inertia 
Impulsive load Wx=M.g=2KN, Wy=-20KN over 2.77msec; 
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Example LOX seal, M=50kg, 26kcpm, DP=37.5MPa, C*=129um T=110K 
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Figure 1 5. Transient LOX seal eccentricity due to step and ramp 
loads. Static load Wx=1 0,000 N, Wy=0 N. 
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Example LOX seal, M=50kg, 26kcpm, DP=37.5MPa, C*=129um T=1 10K 



[N]Aj eojojiees 
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Figure 18. Transient seal force Fy due to impulse load Wy. 

Wx=M.g=490 N, Wy=10 kN (impulse), to=0.57 ms. 
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Figure 19. 1 


Example LOX seal, M=50kg, 26kcpm, DP=37.5MPa, C*=129um T=110K 
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Figure 20. Transient bearing force due to periodic load 
Load 5KN synchronous (433.3 Hz) 
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Load Wx=5 (t/T*) KN, T*=2 msec, 0.2 msec, Wy=0 N. 









X [m] 

Figure 25. Transient journal locus due to impulsive load 
Wx=M.g=490N, Wy=5kN exp(-t/to), to=0.5ms 
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Figure 27. Tn 
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time (secs) 

Figure 30. Initial X-forces for unbalance U=40um, Funb=19.71kN 
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Appendix A. Numerical Evaluation of Unsteady Bulk-Flow in Fluid 
Film Bearings. 

Consider the unsteady flow in the thin film lands of a fluid film 
bearing and governed by the following bulk-flow equations (Yang, 1992): 
Continuity: 

-gl-lp/m,) * ^(pm . 0 <1, 


Momentum: 


Energy: 





( 2 ) 


4 (pHr)+ 4; (p ^ r) 


+ j* 

H 


U 2 X +U 2 y + 


+Q, = 

QRU, 




dp 

i- as 7 


+k 


l&R 2 


J \ 


-U r QR\ 


+Qr” 

2 dX 


(3) 


where i,j=x,y. circumferential and axial coordinates. 

X x =nR; Ay=0, journal surface speeds. 

kj,k,, turbulent shear parameters; and fluid parameters 
C p =C p ( P , T ) ; Pi = Pi( P # T ) , 

P = P ( Pf T ) ; p=jj(P,T) 

with 0,=^ (T-T B ) +hj(T-Tj) as the heat flow through the bearing and 

journal surfaces. 


The film thickness is typically given as: 

H = c(x,y) + X(t)cosQ + y'{t)sin0 

where X(t) and Y(t)are time dependent journal center displacements 

Dimensionless equations of motion are given as: 

Continuity: 

+ (Phuy) + O-^-(P^) =0 


Circumferential-momentum: 

, - - , 

~ h ^c = +Res ~§x {phUx) +i?6p 15 x {phu * 2) (p-hUjpUy) 


Axial-momentum: 


(4) 


(5) 


(6. a) 


-h|P = 

dy 


(fc y U y ) H 




-l^Uy) 


Re \~k (pi3U * U X } + ^ (P hu y 2 ) 


<6.b) 


where , 
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P=P/p.; 

u x=U*/V.; 

c=WR/V.; 

h=H/c. 


u y=U y /V,; 

P=(P-Pa)/Psa 


X=X/R; 

A=fiR/v, 


y=Y/R; 

V„=c 2 Psa//j.R 


t=r/o) 


“ rotational 


MS3mtoiL ^ iaU ^^ 

control voS"‘“cv"K“i'^“/ u “ <«/£> ls integrated on the axial velocity 


ih^y . ff g% £ l7h 

Oil 

+ l[ Re;l ^ Cphu y 2 ) \dxdy 


( 7 ) 


consider the following approximations 


for the different terms i 

O 

lf~ b ^y <ilaiy * f-bp [t>P] n B dx - 


assume the nrpflanro ^ 

uniform film thickness evaluateS^th ° V6r the 8 a "d n faces and 

evaluated at the center of the V-cV r«V d av ®rage 


h p (P a -p P ) 6x v 
fa 

cv. Let, 


m equation ( 7 ). 

(8. a) 




-e., assume average film thickness, viscosity 


And, 


(8.b) 

and shear coefficient k, on V-cV. 

9 il 

(I'd* ( '^ bu - u y ) dxdy = f(P hu x u y )* dy . fit 

w’faSrS St f v-cv! ir ^ Umferential 


Also, 


(9hU x U y) *lyr (8 c) 

flow ( P hu,6y*) across fche east (e) ^ ^ 


9 n 

[fsp <<‘ h “yUy)dxdy . }tph Uy U,)° dx - (fh Uy U y) y ix 

of S iU r v^. SS “ U1 « nniform acroa, the 


(8 .d) 

south (s) and north (n) faces 
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UVPTcvs 



Circumferential velocity 

control volume Axlal veloclt Y control volume 

Figure A.1. Control volumes for integration of equations 
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(8.e) 


i/a? (?hu r } = ^// (phu y ) dxdy 

w a 

since the control volume size is fixed in space. 

Thu., eguation (7, can be written over the control vol„„e , S! 

VlP.-Pe) Sx' . {^]u„ If 6y v . )<* d y 

+ [(P Au xU y ): t>y v +(phu y u y )" e Ixjte; 

integration of the continuity eguation , 1 , on the v-cV give.: 

(phuj; t,y . (e toy); a*- . ,^}f { j h)dxdy . 0 

far j* 


^ the fl °" " teS aCt0 » 3 the volume face, be expr.e.ed 


F i = { P hu J 9 *y v ; K = (phuj - 6 y v. 

F * = ( P hu y )n *> xV -' K = (pAu r )* 8x-; 


Then equation (10) is rewritten as 

F v + F„ - F s + o-^-fJ (ph) dxdy = o 

V 5 


( 9 ) 


( 10 ) 


( 11 ) 


( 12 ) 


?h h i C 5- e cv! abliSheS 3 bulk - f l- balance including the rate of mass accumulation on 

of Lai^ic^r ^^^"Leschzir^r* 8 (1*9 This^ ^ sch^ treated USing the “P^d scheme 

velocities based on whether ihefl a ^-tlon of 

For example 

,0hu.u,r 5 y . f.»;. | p. V, if <0 where P.. p^.-ay 
and can be written as 

F * U ° = lF ” ° lv * " l ~ F *< where la, 0 J=max (a, 0) 


Thus , 


it can be stated that the momentum fluxes 


are equal to: 
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( 13 ) 


(p*u x u y ) " 6y - = F v u? = lF vl olv w - l- Fv , 0}v p 
(phu x u y ) • fiy v = F 9 u y e = I F„ OlVp - l-F„ Olv g 
(phu y u y ) • bx* = F s u y = lF„ 0]V S - l-F s , Qlv p 
(phu y u y )" fix- = F„Uy = [f„ ( Olv, - l-F nl 0K 


The momentum flux difference 


terms in equation 


( 9 ) 


are written as 


Let, 


(phu x u y )Z&yv 


(phu y u y ) n a 


5xV =F » U y ~ F wUy + F n U? - F a Xly 

= V P (lF »' 0U-F v , oMf„, oH-F s , ol) 

- I-F e , 0}v e ~If w , Ol V w -l-F n , Olv w -lF ff/ 0lv s 


a/ = Re p I-F„ ol ; 
a » = Re p I ~F n , Ol ; 


a " = Re p IF v , Ol ; 
a s = R e; lF g , Ol 


and using the identities: 

Ia ' 0l= 1 ; l-a.0l- + 

the following identity is obtained 


F * +l ~ F e> Oj-F v HF v/ 0 ll+F n +l-F n , ObF a +l-F g , 


«' 0JI “ ( F ,- R v +R „-F a ) 

[ a/+a*+a/ +a w v ] 

KOp 


And using the continuity equation (12), 
F »^~ F e> ol+F n +I-F n , o1-f 9 +[f 5> 0l=-o 


-fcff (ph)dxdy 


+ 



where nb refers to the 
Substitution of (16.b) 


neighbor nodes (e,w,n,s). 

into the axial momentum equation (9) gives 


h P v (p a -p P ) fix - (^)*x*y*+Re.£ff(p bUy)dx fy 

+R e;[v p (If,, 0M-f v , oHf„, oH-f s , ol)] 
~ R e; a-F„ Ol V B -lF vl Oiv w -l-F n , Olv^-f f b , Olv s ) 


and substituting equation (16. b). 


(14) 


(15) 


(16 .a) 


(16. b) 


( 17 ) 
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(18) 




VpRe *~&tf. f dxdy 


Hence, the difference form of the momentum equation is obtained as 


h P t Ps-Pp ) V+ g a ii V nb = 




Vp 


( ph) £0x v fiy V ]-A v p 


( 19 ) 


implicit Tchem b e ^s^sey fof the unsteady term, d *^™ tive has to be found. An 


dv P _ Vp-v; 

~ST “ ~ST~ <*> 


where Vp - v P is the axial velocity at time t+At, and 

V P = V P is the axial velocity at time t, respectively. 

evaluated at* time't+At? 3 |vel ° cit '' *"« Pressure, i„ equattou ,19, are to be 

Finally, the discrete form of the axial momentum transport equation is given as, 


p/ (p,-p P ) 6x v +g a^vj + 


Re a (ph) v P t>x v by v 
^ v P = a/ v p 


( 20 ) 


where : 


JM] 

\ h ) l 


bx' 


t>y v+ + Re a (ph) v P 


nh)r a * v 


At 


I ntegration of Circumferent ial Momentum EmiaHnn 

control^olume^fU-cV)' shown^TrTf igure^A^i m °and t us^ r ^ n .h 0rt ‘ 6 '*’ °" ‘he 

simplify some terms (same as for the ' v-cV ^ h ®, co " t:Lnuit y equation to 
algebraic equation: V e( I uatlon ) leads to the following 


where: 


hp lp P -p E ) 8y u + a2u g +aXu w *a2u s +a£u a 


kjV A 


u: = 


h 2 
U at time t 


bx u 6 y u +4^ (0/3)2 


£7 (ph)2 &x u by u u P = a P u a P 


( 21 ) 
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3p - E a ^> + (nr) p 5xU fi r u+ Re s (ph) p 

a e ~ Re p l-F e , Ol ; a* = Re p I F w , oj ; 

a " = Re p I -J a/ Ol ; a s u = /?e p ’ If„, o| 


= (phu x ) 3 by u ; 
F n u = (phu y ) n bx u ; 


•f’r = (phu x ) v by u ; 
F s U = (phu y ) 3 bx u ; 


The difference equations for momentum transport are summarized as 
Circumferential-momentum (on U cV) ; 

hp (Pp~P E ) by U+ ]P a nb +S^ u Up = a P V p 

Axial -momentum (on V cV) : 

hp (Ps-Pe) V^+sZ+S^v; = a/Vp 

where : 


a P = 


h 


bx v by V +Y +5/ ; S t v = Fe s (p/a) V P jZ 

At 




u fiy u + r + S T “ ; 5 X U = i?e s (p/a) “ 


At 


A 

h 2 


u 


At 


s; = o 


in general: 


with 


r= u or v 


«/ = ^ e p I-y/, ol ; 
ay = J? e ; [-F/, ol ; 


•F 1 / = (p/aii x )£ 6y r ; 
F / = (phu y )* bx z } 


< 3 * = /?e p |f', ol ; 

a / = /?e p [f/, ol 


= (P^ a *)* 5y r ; 
= (P^Uy)^ bX T } 


( 22 ) 


(23) 
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CV) Shown i/figSreTl^eads l^ ation < 5 > ° n the pressure control volume(P- 


where the flow rates 


F * F v +F n~Ff +o J J ( ph) dxdy = o 

W 8 

across the faces of the p-cV are: 

F » = (phu x )%&yP . (ph)*Up&yP ; 
F " = (p hu x )ftyP . (p h)*u,fiyP ; 
F ” = (phu y )ft x p . (p A )J^ x p . 

F ' = <P>*u y )5»x^ - (phjfVpSx^ 1 


(24) 


(25) 


The unsteady integral tern, is approximately given by 

^//(Ph)dxcy^ <pi,,J Jyp 5rP . [p5 ^ * A p |?J 8jt „ 5yp 

Which implies that the size of the p cW • 

° f h a " d f “ a “ if »™ withlft^*cL?roi S vlS. ln SP *° e ' a " d that vein., 

T e f£il b llH control ™-ee the f lo„ balance on 


F ?-FS+ F Z- F ?*obxP iyp (pP + hP 3p?l 

l * p ~W\ 


= 0 


(26) 


Since h 


e (x,y) 


+ x(r)cos0 + y ( t) sin0, then 


dhp 

dr 


COS0. 


^(T) 

0T 


sin0 o 


where the time derivatives of x an h - 

bearxng equations of motion. Also? y *** 9iV6n by the ^^tion of the rotor- 
_ Vp-Vp h 

~3T " — Si~ where P=P(t)and p=p ( t+A T ). 

The continuity equation in the P C V is then written as 

F ?-F?+ F n P - F ?+06 X P 6yP |p P + =0 (27) 


using the pressure correction method, .launder and he.chriner, 

U * U + U 1 ; V r =V r +v r/ ; p = ~ ^ , 

u ' . Vy ***"?* f are”cor roc t ion ?Sids 1C " S < 22 a " d 23 >' «■» 


(28) 
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Substituting ( 28 ) 


into the momentum 
h p (P P -p B ) 6 y u + 

(Pp-Ps) fiy u + 


equations (20, and ( 21 ) leads 
Unit + Sp + S x U Up 

^ a nb u nb ~ ZpVp + a/iip 


to 


then. 


V (p;-pi)6y“ + g aj u U 4 = ap u u 


and i 


identically 


Let, using the SIMPLEC 


procedure (Van Doormaal & Raithby, 1984) 

E a *> “4 “ £ a£ uj, 

E a 4 v / 


Equations (30) then become: 


where. 


Let, 


u'-*L 

(Pp~Ps) 


a p 


; dp - hp6y 


(Ps~Pp) 

* dp - hp&x 1 

a; 

~E a 4 


S,», 


Re a (ph) u p A* u fiy u 


At 




, 4 ' r «y r *£ai*^ ; ^ = l?e s (ph)-^16y 

At 


D P U = 


^>4 

d P v 


D P V = 


e «“°“°” t (2”) 0 |i™' ! COrreCti0n field9 -®1u.tio„ (28), i„ t<) th . contlr 


(29) 


(30. a) 


(30. b) 


(31) 


(32) 


173 



where 

F e = (phujy)* = (pjj)* u p b y p 

F s = ( phu y bx) p s = Cph) p v p 6xf , etc 

and 

F * = (pWeUpty* = Cph)iD P u (p^)^ y p 

a s (Pp“pi) ; a/= (“ph)^6yP p p u 

^ = ^ h )> f p*x p - <p«ftp r (p;- P ;)axi> 

= a /(pi“Pp) ; a/= P p v 

_ (p v -pp) 

P/j = a>f (pp-ph) 

Let 

Ff-Ff+Ff-Ff+o&xe 5 y p jp^ p+/lp j»|j5j , j. 

So equation (33) becomes 


(34) 


(35) 


(36) 


F'-fI+f'-f' = -s; 


(37. a) 


or 


or 


a E (Pp P B ) ~ag (Pu~Pp) +ag (p'p-ph) -ai(p' s -p' p ) = 


a * P ' p = ^ a &> PLb ~ s ; 


-s; 


(37. b) 


(37. c) 


where 


a/ - (a B +a„+a N +a s )P = Ya^ 


(38) 


the momentum eolations' are aTscfsatisf ^° ntinuit y is satisfied. Then 
achieved. fied and conversion of the flow field is 

Once the p£ field is obtained, correction to 4 -u 

equations (27). Typically the pressure is unde^laxefal fi6ldS “* d0ne USing 

r->«w _ _okd , cu 0.9 r 


p- = p o« + , 

u”" = u oU + u- 


v“* = V°“ + V' 

With or, as a parameter typically less than one. 


( 39 ) 



D iscretization of Rnerav Tran sport Eguation 


2*«!2TiSgr equati0 '’ f ° r tUrb " le " t ‘“'•“-“O" energy transput is given 


where. 


_ l 'PW * (phu y T) ] ♦ £Skg M 

■ M<ef M|p . 


(40) 


y=Y/R; t=r/u 

U y~^y/U»; A=0R/U.; T=T/T, 

p=(P-Pa)/Psa; h=H/c 


and. 


X=X/R; 

u x = U„/U.; 

0=4) R/U.; 

P« = P, T. 

is the Ecke rd number, 

Re. P =_(?;/^.)« J ^ t hi ^ squeeze 3 f i lm Rfynold3 n n R S n ° ldS number ' 

Qs hetT-TeJ+hjd-T,) is the heat flow to Searing^nd journal surfaces 


Define the following source terms. 




(41) 


= irU + h[ k * {u * +u }+i u A) ) 

in figure A.l leads^o* 311 <40) ° VSr the tem P erafc ure control volume (T-cV) shown 


c TCi E 


Re p — 


I / ( ?hu x T) l dy+ j CphUyT) n s dx I 


+ -jr-QAx&y = (5 1 AxAy) Tp + S 2 AxAy 

implementation of the upwind scheme for the thermal flux terms gives: 

(phu x T) ’Ay = F'T, = I F 9 , Ol T p -l~F 9 , ol T v 


(42) 


(phu x T) v Ay = = If w , ol T w -[-F v , o] T p 


(phu y T) n Ax = lF„, Ol T p -\-F n , 0 J t n 
C phUyT) *Ax = [f b , ol T s -l-F s , ol T p 


where 


F » = (p My) a ; 

F * = (PD*Ay) v ; 


F n = (phUyAx)„ ; 
F„ = (pAUyAx), 
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«3i„ g these sessions the LHS of eguation ,42, is rewritten as, 


or 


LHS *2. = 


= C P Re \ 


where 


!}„l ''““"f * qUf„ oJ]f p 


CpRep ^ __ 

e c °I 7 ’ff + ITF v( o|r w +I-F n , oJr v +ILF o , ofr s ] 

+ Fp^ ep _ \F -p + p -PYT 4. ^ e P m * 

E c i 9 v F o F s ] T p + —^—QgAxAy 


LHS, 


= C P Re - 


'c a IJJ^ {phT)dxdy + £>*>) T P 


a nbTnb 


+ — i F ep rp _p + „ Re P 

E. r* *»* F a F s] T P + ~B-Q b 


AxAy 


(44) 


• L§^l-F.. cl ; a. . 


ol ; 


ol , er-^^,0, 

The diecrete for. of the continuity eguetion on the T CV gives 

= -°f f (ph) dxdy 

TCi 01 

Substitution of equation (46) into (44) gives 


LHS, 


42 


= < z£2i 


tff/J - sg 5 >// £ <ph, drdy 

rev 

+ d2 a nb) Tp - Y a nb T nb + ~-Q a AxAy 


(45) 


(46) 


" ■ S K ! fL l>h % dxdy * E a ») T v - £*»?*, * 

F c f^®// ~§r~ dxdy ~ Re P a f f^ (ph) dxdy^ 


Q s AxAy 


(47) 


^ e> ^ e P a # the last two term on the RHS n<F 1 ~ 
add to zero. Then, he ms of the previous expression is 
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LHS i3 = Sp? e s f r-h dT 


Ef!h h ^ dy * - £a^ ♦ 


OgAxAy (48) 


The integral for* of the energy transport equation (42, is now expressed as: 

^rfP h ^ dxdy * Q>») n 


-Q g AxAy 


(Sj^AxAy) T p + S 2 AxAy 


(49) 


Let, 


— = r f~ t p 

5t At 


Tpihg+hj) - ( h B T B +hjTj ) 


form of equation (49 ^becomes 1 ; 6 ** th “ previous time step. Then, the discrete 

^ 3nb+ ~ P E C ~TT AxA y * ~ (hg+hj) AxAyJr p - (S 1 AxAy)T p 
- ^ 5 2 AxAy + ^f^fAxAyr; + ^(F^F^AxA y 


(50) 


The algebraic energy transport ^theory ^equation is finally written 

3pC Tp = a * T * + a " T * + «» T w + a/ r* + 5/ 

where 


as 


V - -SSit-F.', 0] ; a> « . e /^| r . 

_ C E c 

at - sgzit-Ft, oi , ,« . 


I^V , oi ; 


a/ = 0] 


(51) 


(52. a) 


a P - a/+a h f+a w tr +a/+5^+l5p^ ( oJ+5p^ 


(52. b) 


Sc = 5a+F C2 +S C j+S c <+I-5 P2( 0] T‘ 


(52. c) 
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= 


Re' - _ 

~E~ A B +hj) AxA y 


t n _ h J 


* = c p Re s pA, 


?P3 


~AxA y 


s " 4jtf4 
<? t __ A , . . 4 /Jp 

“ 2* A p(Pa-Pj Ay 


O ^ , T~“ — 

° - ~E~ A B T B+hjTj ) A XpA y p 


t = Pj A, 


X lr AxAyT r 


tern,, evolving it correspond to unsteady flow conditio, 
' ° thS B ° UrCe fe RHS if ,- Sl4xly)<0 . 


(52. d) 


(52. e) 
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INTRODUCTION 


The classical lubrication theory, derived originally by 
Reynolds (1886), assumes that the lubricant is a single-phase, 
isoviscous, newtonian fluid. However, this assumption is commonly 
violated in practical cases, in some processes like metal forming, 


a mixture of water and oil is used as lubricant and provide an 
additional cooling mechanism (Al-sharif, 1992 ). Dust and other 
impurities difficult or expensive to avoid in practical 
applications will turn the lubricant into a solid-liquid two-phase 
mixture (Chamniprasart, 1992). In devices like squeeze film dampers 
entraintment of air, or air dissolved in the oil, will convert the 
lubricant into a gas-liquid two-phase mixture (Zeidan and Vance, 
1989,1990). Operation close to the fluid critical point (like in 
some damper seals using cryogenic fluids) may result in fluid 
vaporization and the consequent formation of a single component 
two-phase fluid (Yang et al., 1993a). In any case, the 
transformation from a single-phase to a two-phase fluid implies 
considerable variations in the overall material properties of the 
fluid such as its viscosity and density. Thus, a pronounced change 
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in the performance characteristic of the fluid film bearing element 
is expected under two-phase flow conditions. 


LITERATURE REVIEW 

The study the two-phase flow phenomena in fluid film bearings 
such as seal, dampers and bearings commonly used in turbomachinery 
is of current interest. Therefore, this brief literature review 
focusses on the behavior of a mixture of a liquid and a dispersed 
compressible fluid (gas or vapor) . A two-phase flow condition is 
attained when the critical point of the gaseous (compressible) 
phase is reached. A gaseous phase may be present in hydrodynamic 
bssrings and dampers by two processes: fluid vaporization (vapor 
cavitation) and gas trapped in the fluid (gaseous cavitation) . 

There have been two main streams in the study of two-phase 
flows as a direct result of their practical applications. The first 
of them considers the fluid as a two-component mixture with no mass 
transfer between the constituents (for lubrication with emulsions 
or for contaminated lubricants) . The second trend of research 
regards the fluid as a one— component , two— phase fluid where the 
vaporization of the liquid requires mass transfer between the 
phases (face seals and annular seals) . The following is a review of 
the relevant research on two— phase flow under these two approaches. 
The continuum theory of mixtures 

One of the approaches used to model a two— phase fluid is 
through the continuum theory of mixtures where the fluid is 
regarded as a homogeneous mixture of two immiscible components. 
This method has been applied to lubrication problems in order to 


180 



analyze the behavior of oil on water (0/W) emulsions used in metal 
forming processes. Sahara et al. <1988, present an experimental 
study on the droplet behavior on 0 ,W emulsion lubrication on a 

"ta 1 r ° lling observations show that the oil 

droplets which penetrate into the elastn h V H 

elasto-hydrodynamic-lubrication 

zone (EHL) form an "oil pool” i e an oi i 

f / i.e. an oil zone containing water 

droplets. In other words, a phase inversion occurs producing a 
water on oil (W/ 0 , emulsion with a higher mixture viscosity in the 
KHb zone. The extension of this oil pool was found to be influenced 
by the on concentration and the roller speed. However, at high 
roller speeds, coarse oil droplets cannot enter to the EHL zone 
because of back flow in the water phase. Only the fine droplets 

moving surfaces enter the EHL zone producing a very 
fine 0/W emulsion. 

Steri and co-workers (Al-Sharif 1992 , Al-sharif, et al 1993 
Wang, et al. 1993 , apply the continuum theory of mixtures to the 
lubrication problem with liguid-liguid emulsions. l„ A l-sharif<s 
work, the continuum theory of mixtures developed by Atkin et al. 
(1976) and Muller (1968) are thoroughly discussed. The main 
assumption of this theory relates to the ..eguipresence principle" 
Which states that every spatial position of a N component mixture 
IS occupied by N particles, one per each component. Conservation 
for each component are derived and coupled by the mutual 
interaction of the components. The authors develop a set of 
constitutive eguations, which when combined with material functions 
for viscosity, thermal conductivity and Helmholtz free energy lead 
to the mixture energy transport equation. The lubrication 




o , 
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approbation (thin films, is applisd to the mass and 
equations render an extended Reynolds equation for mixtures. 
When applied to a journal bearing, the model shows some results 
that agree qualitatively with Nakahara et al. (1988 , observations, 

e*g. a "filtration" phenomenon occurred in the mnu 

iJ7ea ln the converging region of 

fluid film. Due to the difference in viscosity of the 

components, the one with higher viscosity flows easier to the 

minimum film thickness section of the fluid , 41 • 

fluid film increasing the 

mixture viscosity. The model also predicts thar 

u predicts that the maximum peak 

earing pressure increases as the oil content increases. 

in Wang et al. (1993), the continuum mixture theory is applied 
to elasto-hydrodynamic conjunctions. The model intends to represent 
emulsions with suspended deformable particles small enough within 
the minimum film thickness. hn expression for the film thickness as 
function of velocity, load and material properties is obtained. 
Results of the proposed model applied to a oil-i„-„ a ter emulsion 
show that the film thickness varies linearly with the inlet volume 
fraction of the oil. An increase in the concentration of oil is 

observed close to the minimum film thickness , • , 

lckness (oil pooling) . Under 

this phenomenon the film thickness in the EHD • .. 

tne ehd conjunction is 

defined by the viscosity of the oil. 

The next advancement in the mixture theory approach is related 

to gas-liquid mixtures. This type of mixture usually termed as a 

bubbly” mixture may be oresent in =, • 

V be present in applications such as fluid film 

bearings like seals, squeeze film a.m 

9 6 fl1 " 1 dam P e rs and journal bearings 

vented to ambient conditions . 

conditions and under dynamic loading. Tonder 

(1975, 1976a) presents a model for a 

a homogeneous mixture of bubbles 
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and lubricant where the size of the bubbles is small compared to 
the film thickness and with a low bubble solubility. The model 
neglects temperature changes and viscosity variations due to the 
presence of bubbles, i.e. it regards the mixture as isoviscous. 
This model shows that infinitely wide parallel-surface hydrodynamic 
bearings cannot carry load when lubricated by a bubbly oil in the 
absence of large temperature changes. The author concludes that 
"the presence of bubbles in an oil does not generate load carrying 
capacity, it can only modify one produce otherwise." Tonder also 
shows there is a gain in load capacity in externally pressurized 
bearings when bubbles are trapped in the lubricant (to 10% in gas- 
to-lubricant volume ratio) . This load gain is proportional to 
bubble content, bearing load factor, and inversely proportional to 
pressure drop. Tonder' s model presents only qualitative agreement 
with load capacity experimental data presented by Cameron (1966) 
and with findings of Tonder (1976b) . 

Khalil et al. (1980) extended Tonder 's model to account for 
variation in lubricant properties due to changes in pressure and 
temperature along the flow. The model combines continuity, momentum 
and energy transport equations to obtain expressions for the 
bearing mass flow and temperature and pressure gradients. The 
expressions derived are functions of the variable fluid properties, 
film thickness, bearing geometry, and of course, pressure and 
temperature gradients. Application of the model to externally 
pressurized bearings shows that the pressure distribution and load- 
carrying capacity improve with air bubble content for both rotating 
and non-rotating bearings. 
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Chamniprasart et al. (1992, 1993) extend the research of Szeri 
and co-workers on binary mixtures to bubbly-oil mixtures. The model 
considers the fluid as a mixture of a newtonian fluid and an ideal 
gas (typically air). However, the analysis does not account for 
mass transfer between component and makes the model unsuitable to 
handle oil vapor cavitation. The analysis when applied to a finite 
hydrodynamic journal bearing shows that the concentration of air 
(gas) affects significantly the bearing pressure distribution. The 
bubble size seems not to have a significant effect on the bearing 
pressure, although it does affect the size of the two-phase region. 
Comparison of numerical predictions with the experimental pressure 

measurements of Braun et al. (1981) show moderate qualitative 
agreement . 

In general the models based on the continuum theory of 
mixtures present qualitative agreement with experimental and 
practical observations. Results depend on how close the modeled 
flow process (laminar flow with isothermal, isoviscous, adiabatic 
conditions) resembles the actual flow. Note that these models do 
not account yet for mass transfer between phases, and therefore, 
they are not suitable for two-phase flows where fluid 
vaporization/condensation is present. Furthermore, a number of yet 
unknown experimental parameters are needed for accurate results. It 
is important to note that all the models referred consider the 
fluid flow as laminar. This is not an appropriate assumption in 
applications with cryogenic fluids where the trend to higher 
rotational speeds as well as the tendency to roughen the surface to 
improve dynamic stability lead to turbulent flow conditions. 
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Vapor Cavitation 

et .1. ,1.M, present an early attempt to treat the two- 
Pbase now problem arising from fluid vaporisation as a homogeneous 
mixture in thermodynamic non-equilibrium. The model considers the 

conservation equations Tor each component and a ..constitutive 

equation of condensation rsi~ 

ensation or of evaporation', specifying the mass 

rate of liquid formation or extraction per unit volume of mixture. 

usion of this sink or source term is the only difference 
with the standard Reynolds equation for compressible fluids 

Results from the analysis show that this added term has a 
significant effect on the static and dynamic bearing force 
characteristics. Unfortunately, according to the authors, the 
equation of evaporation and condensation which governs this 
source/sink term -is not known for some flow regimes and only in a 
rudimentary way manner for others . » 

Feng, et. al (1,86, present density and viscosity models for 

two-phase homogeneous fluids used in ^ 

used ln dynamically loaded 

hydrodynamic bearings and squeeze fii™ * 

a squeeze film dampers. Expressions for the 

mixture density as a function of the volume * . 

e volume fraction and the vapor 

pressure are derived for two limiting cases: full-collapse 

(solubility equilibrium, and non-collapse (no solubility over a 

load cycle) . The authors show that when the flu in v 

nen rne fluid vapor pressure is 

low it has little influence in the gas phase and can be safely 

regarded as zero (absolute) . An empirical correlation for viscosity 

of certain oil products known as the Hayward model ( 1961 ) is used 
to represent the mixture viscosity. 

Pinkus (1990) rises the issue of the 

r the differences in the heat 
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transfer process in each fiow region. In the all-Uguid region heat 
generation is high and its temperature is slightly higher than the 
stator wail temperature. Thus, heat is transferred from the fluid 
to the walls. In the two-phase region the temperature is lower than 
that of the walls and the heat generation drops with the increments 
of quality. The mixture then extracts heat from the walls to 
maintain the vaporization. In the all-vapor section the heat 
generation is minimum, which when added to the low temperature at 
the end of the two-phase region cause a heat transfer from the 
walls and the temperature of the vapor increases. 

Tv o-Phase Flow in Annular Pr essur e Seals ,nH Fanfi 

The theoretical treatment of multi-phase flow in annular seal 
applications is only recent. Several analysis have presented 
modifications to the classical lubrication theory to account for 
the multi-phase phenomena. Unfortunately, due to the lack of 
experimental data and the variety of approaches used, most of these 
efforts show very different results. 

Hughes et al. (1978) study the vaporization of a lubricant in 
liquid face seals. The flow through the seal is divided into two 
separate regions, one liquid and one gaseous with separate 
equations of motion for pressure, temperature and mass flow rate in 
each fluid phase. The two flow zones are coupled by the conditions 
at the interphase surface. Hughes et al. further assume that both 
phases are at the same temperature (vaporization temperature) . The 
model predicts that the maximum flow rate is reached when the fluid 
is fully liquid, and when two-phase flow occurs most of the 
pressure drop takes place m the vapor phase. The temperature 
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variations across the film are found not to be significant, with 
the main cause for liquid boiling as due to ’'flashing" occuring as 
the pressure decreases along the seal. However, the theoretical 
predictions for fluid pressure and temperature show evidence of a 
transition region where both liquid and vapor phases coexist. 

Yasuna and Hughes (1990) present a "continuous boiling" model 
as opposed to the previous "discrete boiling" analysis. In the 
novel model the flow along the seal is divided into three regions: 
all-liquid, liquid-gas, and all-gas. The two-phase region is 
considered as an homogeneous mixture of saturated liquid and vapor. 
Results show the extent of the two-phase region to be considerable 
in all cases, even when the flow is nearly isothermal (discrete 
model) . For nearly isothermal and dissipation dominated conditions, 
the quality of the mixture is very low (close to zero), and 
explains why the mixture behaves as a saturated liquid. Results 
also show that energy convection effects are important and often 
dominant with respect to the dissipation effects for the typical 
range of clearances used in laminar flow face seals. The "discrete 
boiling" model is then inadequate is such situation. 

Yasuna and Hughes (1992) present an analysis for two-phase 
face seals including squeeze film effects and thermal transients. 
The model is based on the "discrete boiling" model referred 
earlier. Results for simulations of step perturbations of the film 
thickness are presented. For step increments the response decays as 
the rotor approaches the initial steady state conditions. For large 
enough film thickness step decrements the seal collapses due to the 
onset of an unstable vibration. For small step decrements of film 
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thickness the rotor approaches the steady state asymptotically in 

a stable oscillatory motion. The assumption of discrete boiling and 

laminar flow are appropriate only for small film thicknesses (low 

leakage rate) and limit greatly vasuna et al.'s model. As the film 

thickness increases, thermal convection gains importance and a two- 

phase zone may appear. For even larger film thickness the flow will 

certainly become turbulent, and at this increased velocity 

condition the thermal transport process may be considered as 
adiabatic. 

Beatty and Hughes (1987) present a turbulent flow, steady and 
adiabatic flow model for concentric annular pressure seals. Here 

fluid flow is also divided into three regions: all liquid, 
liquid-vapor, and all vapor. Numerical results for the interstage 
seal of the Space shuttle Main Engine High-Pressure Oxidizer 
Turbopump show that the leakage rate is reduced by clearance 
reduction, increment of rotational speed, lengthening of the seal, 
and vapor production (due to pressure drop, . it is concluded that 
subcooling of the liquid before the seal inlet will reduce vapor 
production, and therefore, it win increase the leakage rate _ 

Beatty and Hughes ( 1990 , present a different approach in a 
model for turbulent flow annular seals based on the stratified flow 
of the boiling liquid and vapor phases. Stratification of the 
Phases is attributed to the centrifugal inertia due to shaft 
rotation. Governing equations for film-averaged liquid and vapor 
properties are presented. The model assumes the streams to be 
adiabatic and moving at different bulk-flow velocities, with the 
vapor layer closer to the shaft. Variations of the leakage rate 
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with speed, seal length, vapor generation, etc, are similar to that 
predicted with the previous "homogeneous mixture" model. However, 
the stratified model predicts larger leakage rates and it is less 
sensitive to rotational speed than the homogeneous— mixture model. 

Hendricks (1987) presents a unique experimental study on 
uniform clearance, nonrotating, cylindrical seals for high 
performance turbomachines like the Space Shuttle Main Engine Fuel 
Pump. This report focusses in the measurement of leakage rates and 
P^®ssure profiles for several fluid conditions, and for concentric 
and fully eccentric seal positions. The experiments demonstrate 
that low back pressures may lead to a two-phase flow condition 
within the seal. A model based on one— dimensional, steady state 
isentropic, two— phase momentum and energy equations provides 
expressions for the mass flux and pressure. These expressions are 
normalized to use the corresponding-states theory (Hendricks et 
al., 1973) to compare results using different fluids. Experimental 
results show two-phase choked flow near the seal exit, denoted by 
a positive slope of the pressure profiles in that region. The work 
of Hendricks is relevant because it presents useful experimental 
data for cryogenic fluids, namely nitrogen and hydrogen. Hydrogen 
test data are quite sensitive to the inlet temperature for 
concentric seal configuration. For eccentric operation both, liquid 
hydrogen and liquid nitrogen, are very sensitive to the inlet 
temperature . 

One of the weaknesses of the analyses for two— phase flow 
models referred so far is the lack of firm theoretical grounds for 
the multi-phase transport phenomena (mass, momentum and energy). 
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Gray, W.G. (1975) presents a general equation for multi-phase 
transport of a fluid material property. This work demonstrates the 
importance of the correct definition of the convective and 
diffusive components of the transport process. Gros D'Aillon et al. 
(1981) present a experimental study of two-phase flow through 
nozzles. The measurements, representing the flow under steep 
pressure gradients in nuclear reactors, show the great influence 
that interface mass and momentum transfer have in critical flows. 
It seems that an appropriate treatment of the interface mass, 
momentum and energy transfer could be the key to a more accurate 
prediction of two-phase flows. 

PROPOSED WORK 

The importance of two-phase flow on the performance of fluid 
film bearings is detailed in the previous section. Operation under 
two-phase flow condition is likely to occur in cryogenic liquid 
annular seals operating close to the critical point due to the 
steep variations of pressure and temperature along the seal (Yang 
et al., 1993a) . An accurate theoretical analysis for such operating 
conditions is yet to be developed. 

The objective of this research is to advance a bulk-flow 
analysis for the prediction of the dynamic response of annular 
pressure liquid seals operating under two-phase flow conditions. 
The starting point of this research will be the analysis of two- 
phase flow of a concentric liquid annular seal working with 
cryogenic fluids. As for the structure of the flow, the Beaty and 
Hughes approach (1987) of dividing the flow inside the seal into 
three sections: all-liquid, vapor-liquid, and all-gas ("continuous 
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boiling" model) will be assumed. 

Since the film thickness is very small compared to the seal 
radius, curvature effects can be neglected and the flow domain can 
be unwrapped into a cartesian coordinate system, as shown in Figure 
1* Due l ar 9® axial velocity typical of a turbulent cryogenic 

annular seal, the heat transfer to the walls is much smaller than 
the heat carried by convection, and therefore the flow can be 
regarded as adiabatic (Yang et al. 1993b) . 

The turbulent bulk-flow of a variable properties fluid inside 
an annular pressure seal is described by the continuity, momentum 
and energy transport eguations given by Yang, (1992) : 

Continuity: 

4; {pm ' i * -ft ,pH) * 0 u> 


Momentum : 


Energy : 


where : 


-H- 


dP 

dx 4 


k i u r k j -j- n 
2 1 

H 




H 




* * a ^( p hoiUj) 




if 3 circumferential and axial coordinates 

Ax = fiR ; Ay=0, journal surface speeds, 
turbulent shear parameters 


( 2 ) 


(3) 
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p=p (P, T) ; p=p, (P, T) , 

c p=C p (P,T); 0,=0 t (P,T) for a variable properties fluid. 

These equations are valid for any single-phase compressible 

fluid, and thus can be applied directly to the all- liquid or the 

all-vapor regions. For the vapor-liquid region the transport 

equations have to be specialized to account for the variation of 

the mixture properties such as density, viscosity, etc. The two- 

phase mixture will be considered as homogeneous, and thus, its 

properties will be expressed as a function of the mixture quality. 

Additionally, if the heat of vaporization is significant relative 

to the dissipation energy, i.e. appreciable amounts of heat are 

taken from the fluid to maintain the vaporization, the mixture 

properties will be a function of the heat transfer process (Pinkus, 
1990) . 

The present research will bring special attention to the 
interphase interaction (mass and energy transfer) expressed by the 
mixture quality. The limited understanding of this phenomenon 

accounts for a great deal of the discrepancies among the different 
two— phase flow models. 

As in previous analyses (San Andres, 199i, and Yang et al., 
19«a,b), equations (1-3) win be subject to a perturbation 
analysis to obtain the dynamic response of the seal. The rotor is 
considered to describe small amplitude translational motions about 
an equilibrium position. For such small motions the flow velocity, 
temperature and pressure along with the fluid properties are 
expressed as the superposition of a zeroth order and a first order 
fields. The zeroth order equations describe the equilibrium 
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condition and their solution will provide static parameters like 
leakage, torque and temperature rise. The first order equations 
describe the perturbed dynamic motion of the rotor, and their 
solution will render the linearized force coefficients due to 
journal lateral motions. 
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